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ABSTRACT 

We present measurements of the dust attenuation of Ha-selected emission-Hne galaxies at z = 0.8 from 
the NewHa narrowband survey. The analysis is based on deep follow-up spectroscopy with Magellan/IMACS, 
which captures the strong rest-frame optical emission lines from [OII]/i3727 to [OIII]/15007. The spectroscopic 
sample used in this analysis consists of 341 confirmed Ha emitters. We place constraints on the AGN fraction 
using diagnostics which can be applied at intermediate redshift. We find that at least 5% of the objects in 
our spectroscopic sample can be classified as AGN and 2% are composite, i.e. powered by a combination of 
star-formation and AGN activity. We measure the dust attenuation for individual objects from the ratios of the 
higher order Balmer lines. The li/3 and Hy pair of lines is detected with > 5 in 55 individual objects and 
the Hy6 and H6 pair is detected in 50 individual objects. We also create stacked spectra to probe the attenuation 
in objects without individual detections. The median attenuation at Ha based on the objects with individually 
detected lines is A(Hci')= 0.9 +1.0 magnitudes, in good agreement with the attenuation found in local samples 
of star-forming galaxies. We find that the z = 0.8 galaxies occupy a similar locus of attenuation as a function 
of magnitude, mass and SFR as a comparison sample drawn from the SDSS DR4. Both the results from the 
individual z - 0.8 galaxies and from the stacked spectra show consistency with the mass - attenuation and 
SFR - attenuation relations found in the local Universe, indicating that these relations are also applicable at 
intermediate redshift. 

Subject headings: ISM: dust, extinction - galaxies: high-redshift - galaxies: evolution - galaxies: star forma- 
tion 



1. INTRODUCTION 

Accurate measurements of the dust attenuation in galaxies 
are critical for the determination of the intrinsic physical prop- 
erties of galaxies, and are thus required for studies of galaxy 
formation and evolution. In particular, dust attenuation af- 
fects key observables involving star formation and remains 
one of the largest sources of uncertainty in tracing the star 
formation rate (SFR) volume density over cosmic time (e.g., 
Heddy & Steidel 2009). 

In the local Universe, a common way of measuring the 
nebular dust attenuation is to compare the observed ratios 
of the Balmer emission lines to their expected values in the 
absence of dust attenuation. The Balmer lines arise from 
the recombination and subsequent cascade transitions to the 
second (n - 2) level of hydrogen. The expected rela- 
tive fluxes of the Balmer lines can be computed theoreti- 
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cally (|Menzel&_Bakeij|1937t iBaker & Menzell|1938l; |Seaton| 
1959; Hummer & Storev 1998; Osterbrock & Ferland 200i 
Dopita & Sutherland 2003). These predicted line ratios are a 
function of the temperature and electron density in the nebu- 
lar regions, however the functional dependence is weak and 
thus the line ratio variations are small (~ a few percent). The 
attenuation caused by dust is wavelength-dependent, and it 
is stronger at shorter wavelengths. This wavele ngth depen- 
dence is described by the extinction law (e.g., lO'DonnelH 
IT994tlCalzetti|[T997h . The comparison between the predicted 
and observed line ratios {Ha/H/3 for example) is then used to 
determine the reddening in the nebular regions via the extinc- 
tion law. 

As a physically-motivated and easily-accessible at- 
tenuation diagnostic, Balmer decrements are widely 
employed in studies of the l ocal Universe (e.g. , 
Seaton T979^, Gallagher & Huntei' 198^; 'Calzetti' T99l 
Jansenetal. 2001; Glazebrook et al. 2003; Moustakas et al 
20061; jMoustakas & Kennicutt 20061; ISuU ivan et al 
2001; iKewlevetaL .200L 2002; i Kewlev & Dopita 
2002t iMoustakas et al.l |2010|; IFinkelst eiiTet al I 120111; 
Wiiesinghe et al.l 1201 lallbl) . The Ha attenuation, hereafter 
A(Ha) of normal, nearby, star-for ming galaxies can range 
from 0.1 to 1.8 magnitudes (e.g., iMoustakas & Kennicutj 
l2006t iLee et all ,2009), and is typicaUy ~ 1 magnitude 
in local galaxies similar to the Milky Way. Using star- 
forming galax ies from the Sloan Digital Sky Survey (SDSS), 
iHopkins et al.l (|2003) determine a median A(Hq;) of 1.2 
magnitudes based on the Ha/Hp decrement. Similarly, the 
SDSS study of Brinchmann et al. (2004) finds a mean A(Hci') 
of 1.0 magnitudes with a spread of 0.7 magnitudes for a 
SFR- weighted distribution and a mean of 1.3 magnitudes 
and a similar spread for a stellar mass weighted distribution. 
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B aimer decrements have also been use d to calibrate a num- 
ber of local em pirica l relations (e.g.. iHopkins et aTl 120011: 
iMoustakas et alj|2006l: iKennicutt et al J 120091) . which can be 
used to estimate attenuation in the absence of more direct 
measurements . Such locally-calibrated relations have been 
applied to galaxy samples to z ~ 2, under the assumption 
that the relations stil l hold at higher redshift (e.g., iDale et al] 
IMotlLvet ani20Tlh . 

Clearly, however, the relation between attenuation and lu- 
minosity or mass may change over cosmic time. This expec- 
tation is based on the observation that the amount of dust in 
galaxies has been found to correlate with metal abundanc e 
dHeckman et aLlfigol iBoissier et al]|2004t lAsariet allllOOTh . 
The proposed explanation is that more metal-rich galaxies 
have higher dust to gas ratios and therefore higher extinc- 
tion. The interrelatedness of SFR, ste llar mass, me t alhcity 
and attenuation was also explored by Garn & Best! (1201 Oh . 
who show that the attenuation is most tightly correlated with 
stellar mass (more massive galaxies have built larger dust 
reservoirs), and the metallicity-attenuation relation is deriva- 
tive through the mass-metallicity relation. Yet, the mass- 
metallicity relation has been shown to evolve with redshift 
JSavaglio et al. 2005; Erb et al. 2006; Maiohno et al. 2008; 
iMannucci et al.i 2009. 2010) such that galaxies of a given stel- 
lar mass show lower metallicities at earlier times. Therefore, 
the mass-attenuation and SFR-attenuation relations may also 
evolve in a similar way, and it is important to examine these 
relations with higher redshift datasets. 

For a direct comparison between low and high redshift we 
would ideally want to determine the attenuation consistently 
through the B aimer decrements. Until now, this has not been 
possible for a large sample. Studies of dust attenuation based 
on B aimer decrement measurements have largely been lim- 
ited to galaxies in the local universe, primarily because at 
z > 0.4 Ha shifts into the near infrared, making measurements 
of the Ha/Hp ratio challenging. Alternatively, the higher or- 
der Balmer Hnes HyS (^4861), Hy (A4340), H6 (^4101), which 
are still observable in the optical up to z = 1, may be exam- 
ined. However, observations of these lines are time-intensive 
since they are relatively weak. 

The few studies that examine Balmer decrements in inter- 
mediate redshift galaxies have mainly focused on samples 
which may not be representative of the overall star-forming 
population. For example, some studies focus on luminous 
and ultra-luminous infrared g alaxies (LIRGs and ULIRGs, re- 
spectively) at 0.4 < z < 1.0 (IHammer et al.ll200UrFlores et al.) 
|2004; Liang et al. 20Q3), and thus their findings of elevated 
nebular attenuation (Ay = 1.5, 2.8 and 1.82, respectively) rel- 
ative to typical local g alaxies (Ay = 1.25) may be expected. 
iRodrigues et al.l (l2008h focus on massive galaxies M» > 1 .5 x 
IO'^Mq andfindAv - 1.53 based on the Hy/H/? ratio, in good 
agreement with a complementary estimate based on a compar- 
ison between the infrared and Hf) fluxes (Ay - 1.71). Such 
values are actually comparable to those found locally from the 
SPSS, when a similar range of galaxy masses is considered 
jGarn & Bestll2010l) . Finally, ISavaglio et all |2005) study 28 
galaxies at 0.4 < z < 0.98 from the K-band s elected Gemini 
Deep Deep Survey (GDDS, Abrah am et al.l2004j) . The galax- 
ies have stellar masses mostly between 10^ and 10'" Mq, and 
may be a more representative sample of the star-forming pop- 
ulation during this epoch. However, when the Hy/H/? through 
H8/HyS were measured from a stacked spectrum of the galax- 
ies, the line ratios yielded an attenuation of Ay = 2.13 + 0.32, 
higher than what would be expected locally. 



In this paper, we use the higher order Balmer emission lines 
HyS, Hy, and H6 to constrain the amount of nebular dust atten- 
uation in intermediate redshift star-forming galaxies. Specif- 
ically, we wish to determine whether the average attenuation 
as a function of rest-frame optical luminosity, stellar mass 
and SFR changes significantly between the local universe and 
z - 0.8. This study improves upon previous work by using 
a large sample of Hff-selected galaxies within a narrow slice 
of redshift space down to SFR~ IM© yr"' (i.e., comparable 
to typical SFRs in the local Universe) with deep optical spec- 
troscopy which allows us to detect H/?, Hy and HS. The Ha 
sample selection and spectroscopic observations are described 
in §2. We estimate stellar masses for our z - 0.8 sample us- 
ing spectral energy distribution (SED) fits to the UV, optical 
and NIR photometry. These data-sets and the fitting proce- 
dure are described in §3. In order to make a comparison be- 
tween the attenuation in local and intermediate redshift galax- 
ies, we select a sample of galaxies from the Sloan Digital Sky 
Survey (SDSS) in §4. Our emission-line selected sample is 
likely to contain both objects powered by star formation and 
active galactic nuclei (AGN). In §5 we identify AGN using 
a variety of techniques and place a lower limit on the AGN 
fraction in our sample. The main dust attenuation analysis is 
presented in §6. We present Balmer decrement measurements 
for individual objects at z-0 and z=0.8 as a function of J-band 
magnitude, mass and SFR in §6.1, §6.2 and §6.3. However, 
individual detections are only available for a biased portion 
of the z = 0.8 sample. To remedy that, we create stacked 
spectra from the full spectroscopic sample in bins of magni- 
tude, mass and Hq'-i-[NII]/16584 flux. The Balmer decrements 
measured from the stacked spectra are presented in §6.4. In 
§7 we compare our results with other local and intermediate 
redshift studies, discuss them in the context of galaxy evo- 
lution and examine the limitations and biases in our sample. 
Overall we find that the z = 0.8 galaxies occupy the same lo- 
cus of attenuation as a function of magnitude, mass and SFR. 
The mean nebular attenuation we find is A(Hq')~ 1 magni- 
tude. The stacked spectra also indicate a mean attenuation of 
~ 1 magnitude. The results from the stacked spectra are con- 
sistent with local mass-attenuation and SFR-attenuation rela- 
tions, suggesting that these relations can also be used in the 
intermediate redshift Universe. We summarize our findings 
and discuss future work in §8. 

2. IMACS SPECTROSCOPY 

The galaxy sample used in this analysis consists of 341 Ha 
emitters in a narrow slice of redshift centered at z = 0.8. Deep 
follow-up spectroscopy of these galaxies was obtained with 
the IMACS instrument on the Magellan-I 6.5-m telescope. 
Here we give an overview of the initial near-IR narrow-band 
imaging used to select the sample of Ha emitters, and then de- 
scribe the IMACS follow-up spectroscopic observations and 
their reduction. 

2.1. The NewHa Survey: Narrowband Imaging and 
Emission-Line Galaxy Candidate Sample 

The NewHa Survey is a program whi ch has obtained H g 
selected samples at intermediate redshift (iLee et al.ll201 2a''^. 
Our program has been structured to efficiently capture statis- 
tical samples of both luminous (but rare) and faint emission- 
line galaxies. We do this by combining the near-infrared 
imaging capabilities of NEWFIRM at the KPNO/CTIO 4 
m telescopes (FoV 27^6 x 27 '6) to cover large areas, and 
FourStar at the Magellan 6.5m (FoV 10'9 x 10^9 arcmin) to 
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probe luminosities that are about a factor of three deeper over 
smaller areas. For both cameras, we have designed a pair 
of 1% narrowband filters which fit within high atmospheric 
transmission, low OH airglow windows at 1.18 fim and 2.09 
fj.m. With these two filters we obtain deep Ha-selected galaxy 
samples at z=0.8 (near the beginning of the ten-fold decline 
in the cosmic SFR density) and at z-2.2 (near the peak of the 
cosmic SFR density). 

The work presented in this paper focuses on Ho- emitters at 
z=0.8, which are detected in NEWFIRM narrowband 1.18 fim 
(hereafter NB118) and J imaging of a 0.82 deg" region in the 
Subai-u-XMM Deep Survey (SXDS, iFurusawa et all l2008l) . 
and a 0.24 deg^ region in the SA22 field (a.k.a., SSA22, 
ICowie etal.|[T996l) . Here we summarize the NB118 obser- 
vations, data reduction, and selection method used to produce 
samples of emission-line galaxy candidates that are targeted 
for IMACS spectroscopy. Further details on the NEWFIRM 
NB 118 i magin g co mponent of the s urvey can be found in 
iLv et all (l2Cill and lLee et all (I2012H) . 

Three NEWFIRM pointings in the SXDS (a = 02'' 18"'; 
6 = -05°00') and one in SA22 (a = 22''17'«; 6 = -00°16') 
were obtained in 2007 December, 2008 September and 2008 
October at the KPNO 4m telescope. The cumulative expo- 
sure times for each pointing ranged from 8.47 to 12.67 hr in 
NBl 18 and from 2.40 to 3.97 hr in J. The median seeing dur- 
ing our observations was ~ I'.'l, and varied between I'.'O and 
r.'9, so point sources are adequately sampled by NEWFIRM's 
0'.'4 pixels. Standard near-infrared deep field observing pro- 
cedures and reduction techniques were followed. The 3cr lim- 
iting magnitudes, in apertures containing at least ~80% of the 
flux of a point source, range from 23.7-24.2 AB magnitudes 
in NBl 18 and 23.4-24.1 in J. 

Sources are selected as emission-line galaxy candidates 
if they show a 7-NB118 color excess which is significant 
at the 3cr level and is greater than 0.2 mag. The mini- 
mum of 0.2 mag is based on the scatter in the color excess 
for bright point sources. Coarse corrections for the contin- 
uum slope are applied based on the z' - J c olor, using pub- 
licly available Subaru/Suprime-Cam z' data ( IFurusawa et al.l 
I2OO8) where available. The overall procedure follows gen- 
eral selection techniques commonly used in narrow-band 
surveys (e.g.. iFuiita et al.ll 2003t i Lv et al.ll2007t IShioya et al.l 
l2008t iVillaret al.ll2008n Sobral et al.l 120091) . A total sample 
of 937 emission-line galaxy candidates meeting these cri- 
teria was obtained over the four NEWFIRM pointings. A 
catalog, which contains emission line fluxes and equivalent 
width m easurements of these sources, is provided in iLee et al.i 
(I2012bl) . 

2.2. Spectroscopic Observations 

Deep follow-up spectroscopy of the NewHa NBl 18 excess 
sample was performed in 2008-2009 with the I namori Magel- 
lan Areal Camera and Spectrograph (IMACS; iDressler et al.l 
I2OO6) at the 6.5m Magellan-I telescope. IMACS enables 
multi-object spectroscopy with slit masks over a 27.'4 diam- 
eter area (an excellent match to the field-of-view of NEW- 
FIRM), and has good sensitivity to ~9500A. These two char- 
acteristics make IMACS an ideal instrument for optical spec- 
troscopic follow-up of NewHa NBl 18 excess sources, and in 
particular. Ha emitters at z = 0.80. Our observational set-up 
yields spectral coverage from 6300A to 9600A (correspond- 
ing to rest-frame ~3500 to ~5300A), and captures the strong 
rest-frame optical emission lines from [OIIJ/13727 (observed 



at 6720A) to [OIII]/15007 (observed at 9030A) for galaxies at 
z = 0.8. 

The spectroscopy used in this analysis was obtained over 
three observing runs in 2008 October, 2009 September and 
2009 November, and was mainly focused on sources in the 
SXDS field. Observations in the SA-22 field were carried out 
as time allowed, when the SXDS field was not visible. The 
spectra were taken with the f/2 camera and the Red 300 1/mm 
gratin^3 which has a dispersion of 1.25 A/pix. Masks were 
designed using 7"xl.5" slits. This set-up results in an ob- 
served resolution of 9A, and corresponds to 5 A in the rest- 
frame for the z = 0.8 objects. The WB6300-9500 blocking 
filter was used to limit the length of each spectrum on the de- 
tector and allow more targets to be observed simultaneously. 
An atmospheric dispersion corrector was used in all observa- 
tions. 

Sources were also observed with the same set-up during 
2008 November and 2008 December as mask fillers for a pro- 
gram targeting high redshift Lyo' emitters (PI Masami Ouchi). 
The majority of these observations used a set-up identical to 
ours, but for one mask, a different blocking filer, WB4800- 
7800 was used instead of WB6300-9500. The use of the 
WB4800-7800 filter results in a rest-frame spectral coverage 
of 2670-4200 A range for galaxies at z=0.8. The impact on 
our overall spectroscopic follow-up program is minimal how- 
ever, as this particular mask only included 35 NBl 18 excess 
emitters, eight of which show a redshift of z=0.8. Of these 
eight, five were observed again on other masks, leaving only 
three objects without spectroscopic coverage to the red of 
4200 A. 

Nightly observing procedures were as follows. Bias frames 
were taken before the start of the night. The slit mask was 
aligned on the field using at least 10 bright alignment stars. 
The exposure times for individual science frames were 20 
or 30 minutes. The alignment was checked every hour and 
the mask was re-positioned as needed. Comparison arc lamp 
spectra and flat field exposures were also taken every hour 
throughout the night. The arcs were observed with one set of 
He, Ne, Ar lamps, and flat frames were taken with the quartz 
lamps. At least one spectrophotometric standard was ob- 
served at the beginning of the night, half-way through and/or 
at the end of the night, as time allowed, when conditions were 
photometric. The standard stars were selected to be close in 
air-mass to the science target fields from the ESO cataloJ^. 
The standard stars were observed with a 1.5" long slit mask. 
Arc and flat observations with the long slit mask immediately 
followed. 

Overall, IMACS spectroscopy was obtained for 515 out of 
937 objects in the NEWFIRM 3cr-selected NBl 18 excess cat- 
alog (386/661 in SXDS and 127/276 in SA22). In addition, 
214 SXDS objects which showed color excesses of less than 
3cr significance were also observed to help test the complete- 
ness of the sample. Spectra were taken through 15 different 
slit masks covering 0.62 sq. deg in the SXDS (4 of the masks 
were from the Ouchi et al. high redshift Lyo' program), and 
through 3 masks covering 0.24 sq. deg in SA22. Targets in 
the SXDS were prioritized for inclusion on the masks (via 
the IMACS mask design software) based on their photometric 

There are two 300 1/mm gratings. The blue one has a blazing angle of 
17.5 degrees and a central wavelength of 6700 A. The red one has blazing 
angle of 26.7 degrees and a central wavelength of 8000 A. We used the latter. 



http : / /www . eso . org/sci/observing/tools/standards/spectra/stanlis.h 
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Fig. 1 . — Distribution of tlie cumulative spectroscopic exposure times for 
objects in the sample of 341 Hff-selected galaxies at z = 0.8 which are used 
in this analysis (soHd line histogram). Also shown are the cumulative inte- 
gration times for objects with which were observed during more than one 
observing run and then co-added (dashed line histogram). The medians ai'e 
4.5 and 7.75 hours respectively (dotted vertical lines). 

redshifts dFurusawa et al.l20Tl]) . Highest priority was given to 
sources likely to be intermediate redshift emission-line galax- 
ies (Zphot > 0.7), while low redshift sources {Zphot < 0.7) were 
considered to be low-priority. Slits were placed on 80 to 96% 
of high-priority objects. The remaining high-priority objects 
were not observed due to slit conflicts and generally lie in ar- 
eas with high density of targets. Photometric redshifts for the 
SA22 field were not available to us at the time that the masks 
were designed, and no prioritization was done for the NB118 
excess emitters in that field. 

The cumulative integration times per mask range between 
one and seven hours, with a median of 4.5 hours. Figure [1] 
shows the distribution of cumulative integration times for the 
spectroscopically confirmed Ha emitters that are the focus of 
this analysis. For a sample of 140 NB118 excess sources we 
obtained observations during multiple observing runs. The 
majority of these objects were observed in two runs, how- 
ever 15 of them were observed during three or more runs. 
Most of these observations were incidental duplications ex- 
cept for a sample of 23 z = 0.8 Ha emitters in the SXDS- 
South field, which were intentionally re-targeted to build up 
signal in objects with detected Balmer lines which had initial 
S/N< 5. The cumulative exposure times for the Ha emitters 
which were observed during multiple runs range from 5.5 to 
1 1.5 hours, with a median of 7.75 hours (FigurelTJ. 

2.3. Spectroscopic Reduction 

The IMACS data were reduced using the instrument- 
specific, publicly-available COSMOS software, developed 
by the IMACS inst rument team at Carnegie Observatories 
jOemler et al.ll2008l) . COSMOS relies on an accurate optical 
model of the spectrograph to determine the expected positions 
of the spectral features in the detector plane. The data reduc- 
tion includes the following steps: (1) an approximate wave- 
length solution for each slit is established using three bright, 
isolated arc lines; (2) the wavelength mapping is adjusted us- 
ing all bright arc lines; the mapping is done separately for 
each arc exposure throughout the night; (3) bias frames are 



TABLE 1 
Sample Summary 





SXDS 


SA22 


Total 


3crNB118 emitter 


661 


276 


937 


3cr NBl 18 emitters targeted with IMACS 


386 


127 


515 


Other objects targeted with IMACS" 


214 




214 


Confirmed Ho- emitters with IMACS*" 


299 


42 


341 


Ha emitters with mass estimate 


274 




274 


Objects in spectral stacks 






174 


H/3, Hy detected at > 5(T 






55 


Hj3, H(5 detected at > 5ct 






50 


up, Hy, H(5 detected at > Scr 






34 



^ These are candidate emitters with < 3it significance. 

^ Objects at 0.78 < z < 0.83, independent of their NB118 detection 

significance. 



created by averaging ten bias exposures taken at the beginning 
of each night; (4) flat-fielding frames are created using each 
set of quartz frames throughout the night; (5) bias subtrac- 
tion and flat-fielding are carried out for each science frame; 
(6) sky-subtraction is carried out on all science frames; (7) 
each 2-dimensional science spectrum is extracted to an indi- 
vidual image and its wavelength linearized; (8) the separate 
exposures are co-added. Each science frame is paired with 
the wavelength map and flat field closest in time, typically 
preceding or following the science exposure. The one dimen- 
sional spectra are then extracted using an optimal extraction 
routine ( Horne 1986) which aims to maximize the signal-to- 
noise of the resulting spectrum. An error spectrum is also 
produced. Heliocentric corrections were calculated for each 
mask at the midpoint of the observing night and added to the 
headers. 

Figure 2 illustrates the sky spectrum over the wavelength 
range of our observations, indicating the expected positions 
of the Balmer emission lines over the redshifts range of our 
narrowband selected Ha sample. The figure makes it clear 
that obtaining good sky subtraction is critical for accurate 
line measurement. The COSMOS softw are employs a two- 
dimensional sky subtraction algorithm (Kelson 2003). We 
vary parameters involving the fitting of the sky background 
to minimize the sky residuals. The sky subtraction is done by 
fitting a second order spline to the sky with a 0.75 pixels knot 
spacing. We also exclude the central 8 pixels of the slit from 
the fit to avoid subtraction of the emission lines in the science 
objects. The use of the temporally closest flat-field frame also 
reduced the residuals because it minimized the differences be- 
tween the low-level fringing patterns on the flats and the sci- 
ence frames. Sky residuals were challenging to minimize to 
the red of 8500 A. For this reason, we have visually inspected 
the spectra of all objects included in our work and excluded 
~ 100 spectra (out of the 341) which are dominated by sky 
subtraction errors. 

The standard star spectra were reduced in the same man- 
ner as the science frames. We created a sensitivity func- 
tion for each night using the noao.onedspec. calibrate and 
noao.onedspec.sensfunc routines in IRAFPI. If multiple stan- 
dard stars were observed on a given night, they were all com- 
bined to produce a single, average sensitivity function. The 
sensitivity function was then applied to all science spectra to 

" IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities or Research in Astron- 
omy, Inc., under cooperative agreement with the National Science Founda- 
tion. 
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Fig. 2. — Sky spectram in the wavelength range of our science spectra. The positions of the Balmer lines for redshifts between 0.79 and 0.83, the redshift 
window for which the parent NEWFIRM near-IR narrow-band observations are sensitive to Ha, are shaded light gray. The position of the atmospheric absorption 
A-band is shaded dark gray. 



produce flux-calibrated spectra. This flux calibration, how- 
ever, is not absolute because it does not account for slit losses. 
We side-step this issue by using line ratios in the analysis. 
However, we have estimated the effect of slit losses and off- 
sets in the flux calibration by synthesizing magnitudes from 
the spectra and comparing them to available photometry from 
Subaru/Suprime-Cam in the medium bandpasses NB816 and 
NB921. We find that the spectra are missing ~ 10% + 10% 
of the continuum flux relative to the magnitude measured in a 
2" aperture (for ~1" PSF). 

Redshifts were determined using an automated cross- 
correlation routine based on the SDSS redshift-finding algo- 
rithm (Co ol et al.ll2008h . customized to accept the IMACS 
spectra. The routine uses a minimization to compare each 
observed object to a library of galaxy and QSO templates that 
are linear combinations of eigen-spectra. Each redshift was 
verified by visually examining the 20 best fits as well as 
the sky-subtracted, two-dimensional spectrum to confirm all 
emission features. For 102 objects which were observed dur- 
ing different runs, each epoch of observation yielded a redshift 
measurement. We use these repeat redshift measurements to 
estimate the redshift error of the spectroscopic catalog. The 
mean difference between the two redshift measurements is 
-6.3x 10"'' (-2 km s"') with a standard deviation of 4.9x10"'* 
(146 km s-'). 

Of the 729 objects from the 3cr-selected NBl 18 excess cat- 
alog tai-geted for IMACS spectroscopy, 225 (42) in the SXDS 
(SA22) have redshifts between 0.78 and 0.83, which confirms 
that the narrowband photometric excess is due to Ha emis- 
sion. For 40% of the NBl 18 excess sample in the SXDS 
without spectroscopy. Ha emitters are identified using a color 
classification scheme based on publicly available Suprime- 
Cam RJ'z' photometry, which is calibrated using the galax - 
ies with IMACS redshifts (see Figure 4 in iLv et aLll^Tlh . 
An additional 110 SXDS sources are identified as Ha emit- 
ters in this manner Thus, IMACS spectra are available for 
67% of the SXDS z - 0.8 Ha sample. In the analysis that 
foflows, 74 SXDS sources which have 0.78 < z < 0.83, but 
which show an NBl 18 excess at a significance lower than the 
3cr cut-off are also included. A final color selection to iden- 
tify the Ha emitters among the SA22 NBl 18 excess sources 
without spectroscopy has not yet been performed, but the 42 
confirmed SA22 Ha emitters are included in the analysis as 
well to increase the size of the sample. 



2.4. Spectral Fitting 



In order to obtain continuum- subtracted line flux measure- 
ments, self-consistently corrected for stellar absorption, the 
fluxed spectra were fit with stellar populati on mod els using 
a routine similar to that developed by iTremonti et al.l ([200^ 
for SDSS. This fitting method assumes that any galaxy star 
formation history (SFH) can be approximated by a sum of 
discrete bursts and thus the spectra are fit by linear combina- 
tions of single burst templates. The t emplate library is gener- 
ated from a linear combination of Bruzual & Chariot (2003J) 
models assuming solar metallicity and Chabrier initial mass 
function (IMF, Chabrier 2003). We use templates with ten 
different ages - 10 Gyr, 5 Gyr, 2.5 Gyr, 1.43 Gyr, 0.90 Gyr, 
0.64 Gyr, 0.29 Gyr, 0.10 Gyr, 25.1 Myr, and 5.01 Myr. We 
shift each galaxy spectrum to the rest frame and then construct 
the best fitting model, convolved to the proper spectral reso- 
lution. Internal velocity dispersion of 80 km s"' is assumed. 
The fit inclu des co r rectio n for the foreground extin ction using 
the Schlege l et al.l (119981) extinction map and the lO'DonnelH 
(11994) Milky Way extinction curve. The intrinsic reddening 
is treated as a free parameter using the Chariot & Fall ( 20000 
law. The continuum fit is subtracted from the spectrum, and 
the emission line fluxes and equivalent widths are measured. 
The velocity offset of the emission lines from the systemic 
redshift is assumed to be the same for all lines. The line 
widths are also assumed to be the same. 

In the integrated spectra of all but the most vigorously star 
forming galaxies, the Balmer emission lines typically appear 
to be superimposed on top of the absorption lines of the stellar 
continuum and partially fill them (e.g.. Figure 3). Hence, an 
absorption correction must be applied to recover the correct 
Balmer emission-line flux. The continuum fitting and sub- 
traction procedure provides a self-consistent absorption cor- 
rection to all emission lines. We check whether these cor- 
rections are reasonable, given other values quoted in the lit- 
erature. We measure the equivalent widths of H/3, Hy and 
H6 in the continuum fits and find mean values of 3.1, 2.2 
and 3.4A, respectively. Similar values ha ve been measured 
for local and intermediate redshift samples: Kobul nickv et al.l 
(1999) find that for their sample of local galaxies the ab- 
sorption correction of HyS is between 1 and 6A, with a mean 
of 3+2A:' Mrner & Owenl (12002 ') adopt a morphology depen- 
dent correction assuming that the corrections at H/3, Hy and 
H5 are identical and range from 1.2 to 4.1 (Sa to Irr); and 
ISavaglio et al.l (l2005h determine corrections of 3. 6 A and 3.4A 
for HyS and Hy, respectively, based on a composite spectrum 
of 0.4 < z < 1 .0 galaxies. 

In Figures 3a, 3b and 3c we show three representative spec- 
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Fig. 3a. — A representative spectrum of a z=0.8 Ha-emitter in tlie SXDSS-36609, which has one of the highest H/? fluxes in the sample. This object is also 
identified as composite star-forming and AGN galaxy by the MEx and Spectral type methods (see Section 5). Top: The observed spectrum is shown (black) 
together with the continuum fit (red) as a function of the observed (bottom axis) and the restframe (top axis) wavelengths. Tick marks above the spectrum indicate 
the positions of the Hr5, Hy, and HyS fines. Bottom: A close-up of the Balmer lines US, Hy, and H^S in the observed spectrum with the continuum fit (black and 
red lines) and in the continuum-subtracted spectrum fitted with a Gaussian (black and blue lines). 
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Same as Fi sure [3a| but for the Ha emitter SA22-44387, which has one of the lowest Hfi fluxes. The US emission hne is not detected. 
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tra from our sample. The spectra are ordered by decreasing 
HjS flux. Figure 3a presents an object with a high H/3 flux, 
also identified as a composite SF+AGN object (see Section 5 
), Figure 3b is an object with a mean H/3 flux, and Figure 3c 
is an object with a low Hfi flux. The three objects all have 
the H/3 and Hy lines detected with a S/N> 5 and are included 
in the Hj/H/3 decrement sample considered later in this pa- 
per. While the sky line residuals are visible in the spectra, the 
Balmer lines are clearly detected. The continuum is also well 
detected which allows for a successful continuum fit and reli- 
able subtraction of the stellar absorption in the Balmer lines. 
The distributions of H/3, Hy and H6 line fluxes are presented 
in Figure |4] The H/3 and Hy lines have been simultaneously 
detected in 55 objects with S/N> 5 (left panel). The H/3 and 
H6 lines have been simultaneously detected in 49 objects with 
S/N> 5 (right panel). The line fluxes of the objects in Fig- 
ures 3a, 3b and 3c are indicated above the distributions. We 
note that these fluxes are absorption corrected (as described 
above), but no corrections for slit losses or attenuation have 
been applied. 

3. STELLAR MASSES 

In the local Universe it has been shown that there is a cor- 
relation between stella r mass and d ust attenuati on for star- 
formi ng galaxies (e.g.. lBrinchmann et al. 2004; Garn & BestI 
120101) . Such a correlation may result if star-formation and 
dust-formation are linked such that the build-up of stellar 
mass leads to an increase in dust content and attenuation. 
Using our dataset we can investigate whether any such re- 
lationship holds at z = 0.8. To carry out this analysis, stel- 
lar masses for the spectroscopically observed galaxies must 
be computed. Stellar masses can be reliably obtained with 
spectral energy d istribution (SEP) fits to multi-wavelength 
photometry (e.g., iSaUm et al.il2007t iWalcher et ai]|2011b . A 
wealth of deep multi-wavelength imaging is available for the 
SXDS, and we use these data to derive stellar masses for the 
spectroscopically confirmed Ho- emitters in that field. Specif- 
ically, we use the publi cly available Subaru Suprime-Cam 
SXDS optical dataset ( iFurusawa et al.ll2008l) . combined with 
the NEWFIRM J-band photometry and new GALEX NUV 
observations obtained by the NewHo- team (P.I. S. Salim, GI6- 
005). Comparable data are not available for the SA22 field, so 
we confine the portions of the analysis involving stellar mass 
to the SXDS Ha sample. In this section we describe the NUV 
and optical observations, and t he SED fitting procedur e. The 
J-band imaging is described in iLv et aT] (1201 1) and L ee et al.l 
(l2012b.) . and is also summarized in Section 2.1 above. 

3.L Optical Photometry 

Deep optical imaging in the SXDS over 1.22 deg^ in five 
contiguous fields has been obtained with Suprime-Cam on the 
Subaru 8.2-m teles cope as part of the Su baru Telescope Ob- 
servatory Projects (Furus awa et al.ll2"008h . Imaging has been 
performed in five broadband filters, B, V, Rc, /', and z', to the 
depths of B = 28.4, V = 27.8, Rc = 27.7, /' = 27.7, and 
z' = 26.6, respectively (AB, 3cr, 2"aperture). The three fields 
in the SXDS observed with NEWFIRM are roughly coinci- 
dent with the Suprime-Cam SXDS-S, SXDS-W and SXDS-N 
pointings (Fig. 2, Naka iima et al. 2011) . The reduction of 
these data is described in lFurusawa et aTl (120081) . 

3.2. GALEX Data 

Near-UV (NUV, A - 2316 A) photometry was obtained 
from deep imaging with the Galaxy Evolution Explorer 



(GALEX, [Martinet all 12001 iMorrissev et aH l2007l) . At 
z = 0.8, the NUV band samples the rest-frame far-UV, provid- 
ing strong constraints on the star formation rate and history, 
and additional leverage on the stellar mass-to-light ratio and 
the stellar mass. Our NEWFIRM fields in the SXDS partially 
overlap with three archival GALEX "tiles" (1.2 deg diameter 
circular fields) which have exposure times 28, 26 and 30 ks. 
To increase the overlap and the eff'ective depth, we acquired 
one additional GALEX tile with an exposure time of 46 ks. 
The mean 5 sigma depth for the combined NUV imaging in 
the NEWFIRM fields is 25.3 AB (exposure time -80 ks). This 
is only ~ 1 magnitude s hallower than the deepest NUV tile (the 
EGS field; ISalim et aL 20091 which, as of GALEX Release 
6 ("GR6") has a total integration time of 257 ks. 

At z = 0.8, sources are unresolved in the GALEX imag- 
ing (the GALEX FWHM is 5") and are subject to blending. 
The GALEX pipeline catalogs, based on aperture photometry, 
do not optimally deal with these issues. Therefore we carried 
out a custom PSF source extraction using DAOPHOT. We set 
the DAOPHOT detection and photometry parameters to opti- 
mize the completeness of the detections, and currently gener- 
ate separate flux catalogs for each of the four GALEX tiles. 
Each source list is matched to the NEWFIRM source list using 
a 3" search radius, which corresponds to ~ 3cr positional un- 
certainty. Of the 299 spectroscopically confirmed Ha sources 
at z = 0.8, 85% are detected in at least one of the NUV tiles. 
The photometry from the four source lists is then averaged for 
objects detected in more than one tile, where the average is 
weighted by photometry errors. More details on the GALEX 
data avai lable in the SXDS field is provided bv Salim et 
din prep.h and lLee et al.l (l2012bl) . ISalim et al.1 din prep.h will 
also present a catalog of the de-blended photometry, as mea- 
sured from the stacked SXDS mosaic, and further optimiza- 
tion of the NUV photometry is in progress. 

3.3. SED Fitting 

Stellar masses are estimated by fitting SEDs constructed 
from the multi-wavelength photometry described aboveQ 
The SEDs are fit with a grid of BC03 stellar population 
synthesis models. The methodology is identical to that of 
ISalim et all dlOOTl l2009h and we refer the reader to these pa- 
pers for a detailed discussion. Model libraries are built with a 
wide range of star formation histories and metallicities. Only 
models with formation ages shorter than the age of the uni- 
verse at a given redshift are allow ed. Each mo del is attenuated 
according to the prescription of iCharlot & Fa ll (2000), with 
variable values of both the total optical depth and the fraction 
of the total optical depth due to attenuation of the ambient 
ISM (tv and yU, respectively). The dust attenuation is heav- 
ily constrained by the UV slope which ge ts shallower with 
increasing attenuation (ICalzetti et al.lll994h . However, differ- 
ences in the star formation histories can produce sign ificant 
scatter between the UV slope and dust attenuation (Kon g et al.l 
I2OOI . This is overcome in our modeling because the NIR 
and optical data help to constrain the age. Inter galactic red- 
dening is included via the prescription of Madau e t alj dl996h 
and a Chabrier IMF is assumed. The spectroscopic redshift 
provides the luminosity distance which allows the apparent 
model quantities to be scaled to absolute values. 

This SED fitting is distinct from the fitting of the IMACS spectra dis- 
cussed in Section [23] wliich provided restframe optical (3700 5200A) 

continuum levels for the continuum subtraction. The SED fitting described 
in this section is based on photometry from the NUV to the NIR. 
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Fig. 4. — Distributions of tlie fluxes of the Balmer lines in tlie samples witli individual line detection after the absorption correction. Left: The H/3 (solid line) 
and Hy (dotted line) emission lines are simultaneously detected with S/N> 5 in 55 objects. The vertical lines above the distributions indicate the H/J and Hy 
fluxes of the objects shown in Figures 3a, 3b, and 3c. Right: The Hfi and US (dashed lines) emission lines have been simultaneously detected with S/N> 5 in 48 
objects. The vertical lines above the distribution indicate the Ufi and US fluxes of the objects shown in Figures 3a and 3b. US is not detected in the objects shown 
in Figure 3c. The line fluxes presented here have not been corrected for foreground attenuation. 

The SED fitting performed in this study involves up to six 
flux points (NUV, B, V, i, z, J), their photometric errors and the 
spectroscopic redshifts. For the spectroscopically confirmed 
z = 0.8 Ha emitters, only 20% of the sample is fit with fewer 
than six flux points, and a minimum of four points are used. 
The r'-band photometric points are excluded for the z = 0.8 
sample in order to avoid the re stframe 370 0A discrepancy 
with the stellar synthesis models (ISalim et al.ll2009.) . For each 
galaxy the observed flux points are compared to the model 
flux points and the goodness of fit (x^) determines the proba- 
bility weight of a given model. The average of the probability 
distribution of each parameter (e.g., stellar mass) is the nomi- 
nal estimate of that parameter and its width is used to estimate 
the errors and confidence intervals. 

The final sample of z = 0.8 Ha emitters with determined 
masses is 274 objects. The majority of galaxies in the sample 
are well fit and the median x per degree of freedom of the 
best fitting SED models is one. Seven objects are excluded 
from the sample due to poor fits (i.e., if of the best fitting 
model is >10). Of these seven, one is a broadline AGN and 
one is a possible gravitational lens with two sets of emission 
lines at different redshifts. The remaining five poor SED fits 
are likely due to errors in the photometry. 

Stellar mass uncertainties for the z = 0.8 Ha sample are on 
average 0.10 dex (25%, Icr confidence interval) and are de- 
termined from the widths of the probability distribution func- 
tions for each galaxy. These estimates include uncertainties 
from input photometry and parameter degeneracy (e.g. with 
respect to SFH and dust). Additional systematic uncertain- 
ties may arise from the mod els themselves and t he choice of 
IMF (e g., Maraston 20051 IConrov et al.l 120091; iTavlor et all 
1201 Ih . ISalim et al.. (.20051) show that the masses derived via 
the method used here are in good agreement with those of 
iKauffmann et al.l (12003') and the scatter in the differences be- 
tween the two mass estimates is well matched by the uncer- 
tainties in the two studies. 

The distribution of derived stellar masses for the SXDS 
spectroscopically confirmed z = 0.8 Ha emitters is shown 
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Fig. 5. — Stellai' masses for the NEWFIRM sample of z = 0.8 Ha emitters 
in the SXDS field. The number of galaxies for which stellar masses could be 
derived with uncertainties of ~25% is 274. The masses of the sample range 
from 6 X 10^ Mq to lO'' 
line). 



Mq, and the median is 6 X 10 Mq (vertical dotted 



in Figure|5] The median of the sample is 6 x lO^M©, and 97% 
have masses between 1 x lO^M© and 4 x lO'^M©. There are 
a few galaxies with masses as low as 6 x IO^Mq, and as high 
as 1O"M0. 

4. SDSS COMPARISON SAMPLE 

To Study the evolution of the dust attenuation between red- 
shift 0.8 and today, a comparison sample of star-forming 
galaxies at z ~ is required. Here, we use the data from 
file Sloan Digital Sky Stirvey, D ata Release 4 (SDSS DR4; 
lAdelman-McCarthv et alj [20061) . This dataset has been re- 
processed by a group from the Max Planck Institute for As- 
trophysics and the Johns Hopkins University (hereafter, the 
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MPA-JHU groufPh. The a nalysis pipelines for prev ious 
data releases are des cribed in iBrinchmann et alj ( 120041) and 
iTremonti et alj (2004). Spectroscopic redshifts and line fluxes 
are available for 406,506 objects with 0.04 < z < 0.2. 

To select a local sample appropriate for comparison we ap- 
ply a series of cuts to remove galaxies with low-quality spec- 
tra and uncertain flux measurements. Removal of AGN from 
the sample is discussed in the next section. Sources which 
have uncertain redshift estimates (Z.WARNING> 0) and red- 
shifts at 0.146 < z < 0.148 are excluded. The latter criterion 
is applied because of contamination of the redshifted H/3 line 
from the [OI]/15577 sky emission line. Multiple entries of the 
same object are also excluded using the list of repeats pro- 
vided by the MPA-JHU group. Balmer line ratios are used to 
derive the nebular attenuations of the SDSS galaxies, so we 
also require that Ha and H/J are both detected at a S/N> j^ . 
150,942 objects remain after these selection criteria are ap- 
plied. If the S/N> 5 restriction is also placed on the Hy and 
H6 fluxes, the sample is reduced to 123,950 and 67,261 ob- 
jects, respectively. 

To compare the relationship between dust attenuation and 
stellar mass at z ~ and z = 0.8, the stellar masses for 
the SDSS galaxies derived by Gallazzi et al. (2005, 2006) are 
used. These masses are computed using a Bayesian likeli- 
hood estimate based on fits of the D4000, UJi, H6a + Hja, 
[Mg2Fe] and [MgFe]' measurements. The models account 
for dust attenuation, including negative attenuation due to er- 
rors in the index and magnitude measurements. The stellar 
masses derived via this method for the SDSS sample are in 
good agreement with those of Kauffmann et al] (I2003h . 

To compare the attenuations of the SDSS galaxies and the 
NEWFIRM Ha emitters as a function of the observed J-band 
magnitude for galaxies at z = 0.8, we shift the SDSS SEDs to 
z - 0.8 and synthesize J-band magnitudes using the bandpass 
of the NEWFIRM J-band filter. The SDSS SEDs are con- 
structed usingugriz "Cmodel" tot al magnitude j^y ] and th e 
IDL kcorreci^^ procedures (v4^, iBlanton & RoweisI 120071) . 
which fits restricted SED models to the photometry, using 
BC03 stellar evolution templates. 

Finally, we also examine the dust attenuation as a function 
of the SFR. To calculate the SFRs for the SDSS sample, we 
use Ha fluxes measured from the SDSS spectroscopy. The Ha 
SFRs are not dust-corrected for consistency with the z = 0.8 
sample. The SDSS Ha fluxes are measured within a 3" fiber 
which only samples a portion of the galaxy. Therefore we 
adopt an aperture corrections based on the difference between 
the fiber and Petross ian r'-band magnitudes as suggested by 
iHopkins et alldlOOl : 



2Q-0.4(4^,„-r},,j„) 



(1) 



Of course, it must be kept in mind that the presence of ra- 
dial gradients within galaxies or other spatial variations in 
the SFR may lead to systematic u ncertainties when this cor- 
rection is applied. For example, Brinchmann et al.l (120041) 
have suggested that this aperture correction may be problem- 
atic for galaxies with stellar masses logM, > 10.5 because 

|http : //www ■ mpa- garching . mpg . de/SDSS/DR47| 

The formal line flux uncertainties in the Ha, 1^8, [OIII]/15()07 and 
[Nn]^6584 lines are increased by factors of 2.473, 1.882, 1.566 and 2.039, 
respectively, as recommended on the SDSS data server. 

" http://www. sdss .org/dr3/algorithms/photometry .html#mag_mode 
^'' ihttp:/ /www . mpa- garching . mpg . de/SDSS/DR4/SDSS_inf o . htmlf 
^' http : //howdy . physics . nyu . edu/index . php/Kcorr ect 



they have large bulges which contain little star-formation. 
Discussions of the potential systematics in the estimated to- 
tal SFR resulting from such corrections can be found in 
Kochanek et al. (2000), Baldry et al. (2002), (Gomez et^H 
2003), Nakam ura et all (I2003D . and iPerez -Gonzale z et all 
(.2003) . Using these corrected Ha flux es, the observed SFR 
is calculated using flie lKennicutj ([19981) prescription. 

5. AGN IDENTIFICATION 

The objective of this study is to extend Balmer decre- 
ment measurements of the nebular attenuation of typical star- 
forming galaxies to intermediate redshift. To do this we must 
first identify and exclude galaxies with strong AGN signa- 
tures from the samples. Emission that is dominated by AGN 
light may only represent the small central region of a galaxy 
and not the star forming regions in general. Furthermore, the 
physical conditions of the gas in the central regions of AGN- 
powered galaxies may result in diff'erences in the dust cre- 
ation and destruction processes from those operating in star- 
forming regions. 

For the purposes of our work we are only interested in re- 
moving the AGN with rest-frame optical signatures which 
dictates the choice of methods we use. Our methods can be 
broadly divided in two categories: ones that identify Type 1 
AGN and ones that identify Type 2/LINERs. Type 1 AGN do 
not obey the diagnostic developed to identify Type 2 AGN. 
Thus we first exclude the Type 1 and then proceed with with 
Type 2 identification. 

We consider two methods to identify Type 1 AGN: NUV 
variability and optical line width. Type 1 AGN typically ex- 
hibit broad emission lines which are thought to originate from 
material close to the central black hole. The lines are gravi- 
tationally broadened. In addition, the flux emitted by quasars 
and AGN has been found to vary aperiodically, on time scales 
from hours to decades at all wavelengths. Such variabil- 
ity wa s first reported in connection with reverberation map- 
ping ( Peter sonll 19931: iPaltani & Courvoisieill994t iKaspi et al.l 
1^00, 2007). Evidence suggests that the majority of UV- 
optical variability is driven by changes in the disc accretion 
rates which cause variation in the UV continuum (Welsh et al] 

I26nh . 

In the SDSS sample Type 2 AGN are identified using the 
BPT diagram dBaldwin et al.lll981h. In the NewHg sample, 
we u se the Mex diagram dWeiner et all l2007t lJuneau et al.l 
l20Tlh . 

In the local Universe, a method widely used to dis- 
tinguish star-formation from Type 2 Seyfert and LINER- 
powered emission is the BPT diagnostic, which is based 
on the [OIIIU5C)07/H 6 an d [NIIU6584 / Hq' li ne ratios. 
iKauffmann et al.1 ( 120031) and iKewlev et all ( 1200 Ih have de- 
fined criteria in this parameter space to separate galaxies 
by the dominant source of emission - star-formation, AGN 
(Seyfert 1 or LINER) or composite SF+AGN. We apply 
this diagnostic method to our comparison sample of SDSS 
emission-line galaxies. 

At z = 0.8, however, the Ha and [NII]/16584 lines are red- 
shifted to the NIR and are no longer accessible with optical 
spectroscopy, and thus the BPT diagnostic cannot be applied. 
In its place we use a combination of alternative methods to 
identify AGN based on the overall spectral characteristics, 
blue spectral features, and the physical properties of the galax- 
ies. In the future, it will be valuable to obtain NIR spectra for 
the z = 0.8 Ha emitters to identify AGN with the BPT diag- 
nostic. 
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The MEx diagnostic was proposed to distinguish star- 
formation and A GN emission at z > 0.4 dWeiner et al.ll200'7l : 
I Juneau et al.l 1201 li) . The method relies on the combination 
of [OIII]/l5007/Hj8 line ratio and the stellar mass. The stellar 
mass is used as a substitute for the traditional [N11]/16584/Hq' 
ratio because it is a metallicity indicator for star-forming 
galaxies, and the empirical mass-metalicity relation suggests 
a p hysical connection be tween stellar mass and that line ra- 
tio jTremonti et al.ll2004l) . AGN have higher [N11]^6584/Hq' 
ratios because of the excitation conditions of the gas, but 
they also reside in gala xies with the highest stellar masses 
( iKauffrnann et al.ll2003h . An alternative diagnostic diagram 
which us es the (U - B) co lor instead of stellar mass is pre- 
sented by lYan et alj ( 1201 Ih . This color-excitation diagnostic 
agrees well with the MEx method thus here we have decided 
to only use the latter. 

The objects identified as AGN with the various methods in 
this Section are listed in Table|2l 

NUV variability: The GALEX data, which include six 
epochs of observations spaced by ~year-long timescales, al- 
low us to look for variations in the rest-frame FUV 983 - 
1572 A (observed NUV) emission. One variable Ha emitter, 
whose observed NUV flux varies by 1.2 magnitudes, is iden- 
tified. The optical spectrum of this object is dominated by 
broad emission lines. 

Broad Lines: The routine used to fit the continua in the 
IMACS spectra also provides coarse spectral type classifica- 
tions for the objects based on the best fitting template. Of all 
341 Ha emitters, five are identified as QSOs and/or broad- 
lined AGN. The NUV variable discussed above is also iden- 
tified as an AGN by the fitting routine. The continua of two 
objects identified as hosting broad-lined AGN are not well-fit 
and are excluded from the remainder of the analysis. Finally, 
all of the spectra for the z = 0.8 Ha emitters were also vi- 
sually inspected, and two additional objects with broadened 
lines which were not identified by the routine are found. 

The MEX diagnostic: To apply the MEx diagnostic to our 
sample, we require that both the [Olll]/i5007 and the Hfi lines 
are detected with S/N> 3. Of the 341 emitters, 221 fit this 
requirement. As noted above, a fraction of the sample objects 
do not have mass estimates due to lack of optical photometry 
(the 42 emitters in SA22 and 12 of the emitters in SXDS). 
The objects with poor SED fits (13 objects in SXDS) are also 
excluded. In the end, the method can be applied to 175 of the 
341 z = 0.8 Ha emitters with IMACS spectra. 

Figure|6]shows the MEx diagnostic for the z = 0.8 Ha emit- 
ters and the SDSS DR4 comparison sample. There is an excel- 
lent overlap between the loci of the two samples. In addition 
to this general demarcation of the regio ns predominantly oc- 
cupied by the different types of objects, I Juneau et alj ( 1201 Ih 
have also made public their code which, given line fluxes and 
mass measurements with errors, calculates the probability that 
a given object belongs to one of four groups: star-forming 
galaxy, composite object, LINER or Type 2 Seyfert galaxy. 
We use these probabilities to identify the emission source for 
the objects in our sample (after discarding objects where the 
emission lines are affected by poor sky line subtraction). We 
identify nine AGN and seven composite objects (Figure |6]l. 
Therefore the AGN fraction in the NewHa sample is at least 
5% and at least 4% of the objects are composite. There are 
seven objects which lie close to the line separating the star- 
forming and composite galaxies and have almost equal prob- 
abilities to be part of the two classes. In these cases we select 
the highest probability, which is for a star-forming galaxy for 
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Fig. 6. — The Mex diagram lJuneau et aU2Ql ll) is used to quantify the AGN 
contamination of the sample. The local SDSS sample is split by their position 
on the traditional BPT diagram: star-forming (black), composite (light gray) 
and AGN (dark gray). Triangles show the objects from the NewHa sample 
with lo" eiTors. We require that the [Olll]/i5007 and HyS line be detected with 
S/N> 3. The MEx diagram identifies nine AGN candidates (blue points) out 
of a sample of 1 80 galaxies for which we have both line measurements and a 
mass estimate. We also identify seven candidate AGN-SF composite objects 
(cyan triangles). The remaining objects are normal star forming galaxies (red 
points). 

five of them. 

Stellar masses are not available for the SA-22 sample. We 
use the tight correlation between J magnitude and stellar mass 
to place the SA-22 Ha emitters on the MEx diagram as shown 
in Figure |7] Once the linear relation is removed, the residual 
scatter in the mass is 0.2 dex which we use as an error on the 
mass estimate. In this way we identify three AGN candidates 
in this field. 

The MEx diagnostic is defined at z = and we might 
expect that it will evolve with redshift due to the evolu- 
tion of the physical properties of galaxies. Such evolu- 
tion may be caused by the evolving mass-metallicity relation 
(e.g., Erb et al. 2006), downsizing, additional AGN contribu- 
tion (e.g. Groves et al. 2006; Wright et al. 2010) or cha nging 
physical conditions in the H II regions (e.g., Brinchman n et al.l 
120081: iHainline et ani2009l) . However, further studies of inter- 
mediate redshift galaxies are required to quantify these evolu- 
tionary effects and disentangle them from the possible selec- 
tion effects of current samples. Juneau et al. (201 1) find that 
such eff'ects are likely minor because at 0.3 < z < 1 the MEx- 
selected AGN have X-ray properties consistent with those of 
AGN - they are either X-ray luminous or have a hard X-ray 
spectrum when stacked. Star-forming MEx selected galaxies 
show only a soft signal in X-ray stacks, consistent with purely 
star-forming systems. 

Combining the results from the diff'erent methods yields 17 
AGN candidates and seven composite galaxies (Table |2]i. It 
is unlikely that we have identified all AGN in the sample, 
so this is a strict lower limit on the number of AGN in the 
spectroscopically confirmed z - 0.8 Ha-selected sample. We 
are unable to apply the MEx diagnostic to 47% of our emit- 
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ter sample because we either have no mass measurements or 
the S/N in the lines is lower than our limit. However, we are 
likely to have identified the strongest, unobscured AGN and 
these are the ones which might affect the Balmer line ratios. If 
the 17 objects identified as AGN and the seven objects iden- 
tified as composite are the only ones in the sample, then the 
AGN fraction is at least 5% and the composite fraction is at 
le ast 2%. Th i s AGN fraction is consistent with the findings 
of iGarn et al.l (120101) who also examine the AGN contamina- 
tion in a sample of z = 0.84 Ha-selected emitters similar to 
ours. They identify 5 - 15% of their sample as AGN using 
a broader range of methods which includes NIR color-color 
diagnostics, mid-to-far IR slope and X-ray emission. Again, 
we underscore that for the purposes of this work, we are only 
concerned with removing AGN with rest-frame optical spec- 
tral signatures, which explains why our fraction is at the lower 
end of^ their range. 

We remove AGN from the SPSS sample using the classical 
BPT diagram ([Baldwin et alJ[T98ll) . To apply the diagnostic, 
we require S/N> 3 for [OIII]/15007 and [N1I]^6584 as a min- 
imum requirement for these lines to be reliably detected, in 
addition to the S/N> 5 cut imposed on the Balmer lines in 
Section 4. The cuts on [Olll]/I5007 and [Nll]/16584 remove 
only a total of 2,998 objects, for a final sample of 64,263 ob- 
jects. The most stringent requirement on the sample is the 
H(5 line de tection, thus the sam ple i s effectively H(5-selected . 
Using the Kewlev et alj ( 1200 Ih and iKauffmann et alj ( l2003h 
relations, 6,629 objects (10.3% of the sample) are identified 
as composite, 2,417 objects (3.8 % of the sample) are iden- 
tified as AGN and the remaining 55,217 (85.9 %) - as star- 
forming. In Figure |6] these three classe s are plotted in th e 
parameter space of the MEx diagnostic d Juneau et al.ll201 iT) . 
In order to be able to apply this diagnostic we also require that 
[011]/13727 has S/N> 3 (removing 116 objects) and that the 
objects have mass measurements (removing 17,565 galaxies). 

6. BALMER DECREMENT ATTENUATION MEASUREMENTS 

Comparison of the observed ratio between any two Balmer 
lines to the expected recombination ratio (the Balmer decre- 
ment) yields a measure of the nebular dust attenuation, given 
an extinction curve. A detailed description of how the atten- 
uation is calculated from the observed Balmer ratios is given 
in the Appendix. Here, we assume standard Case B recom- 
bination with temperature T - 10"* K and electron density 
«(, = 10^ cm"^, which results in expected ratios of 2.86, 
0.47 and 0.26 for Ha/^S, Hy/t^S and H(5/H/3, respectively 
(jOsterbrock 1989). We also assume the Calzetti (1997) at- 
tenuation law to convert the color excess to attenuation. 

The present paper is focused on the high order Balmer ra- 
tios Hy/Hj6 and H5/Hj6. Of course, the z - 0.8 sample is Ha- 
selected, and we do have measurements of the Ha line fluxes. 
However, we do not include analysis of the Ha Balmer decre- 
ments in this work because it requires a detailed examination 
of the possible sources of systematic error when combining 
the photometric and spectroscopic data and it is beyond the 
scope of the current paper The challenges in measuring the 
Ha/HjS ratio with the existing data are discussed in more detail 
in Section 7.2. 

6.1. Nebular Dust Attenuation in a Comparison Sample of 
Local Galaxies 

Before examining the results for the z - 0.8 Ha-selected 
galaxies, we review the attenuation properties of the local 
SDSS sample. The Ha attenuation, A{Ha), is computed from 



the Ha/H/3 line ratio, as well as from Hy/HjS and H5/Hj8. In 
Figure[8] the mean weighted A(Ha) is shown as a function of 
J magnitude, stelar mass and SFR. The attenuation measure- 
ments derived from Ha/HyS and from Hy/Hy6 are consistent 
and the trends are in agreement. However, the attenuation 
derived from the HcJ/H/? ratio shows systematic offsets. The 
attenuation based on H5/Hy6 is lower at bright magnitudes, 
high masses and low star-formation rates. Such an offset is 
also reported bv I Groves et al.l (1201 1[) based on measurements 
from the SDSS DR7 dataset and is thought to be caused by 
a systematic error in the continuum fit around the H(5 line. 
The exact source of the discrepancy, however, is still being 
investigated by the SDSS collaboration (J. Brinchmann, pri- 
vate communication). We do not try to correct for this offset 
in our analysis, and in the figures that follow, we will com- 
pare our z - 0.8 sample to the SDSS Hy/Hj6 sample, unless 
otherwise noted. The median attenuations and standard devi- 
ations from the three line ratios are 0.98 + 0.43, 0.97 + 0.82 
and 0.62 + 0.90, respectively. 

The scatter dramatically increases when A(Ha') is based on 
the Hy/Hj6 and H5/Hy6 ratios because of the weakness of the 
lines, the limited wavelength range spanned by them and the 
associated increase in the measurement errors. The large scat- 
ter is one reason for the larger fraction of objects exhibiting 
negative attenuation: while only 0.5% of galaxies have neg- 
ative attenuation when A(Ha) is computed from Ha/HjS the 
fraction increases to 8% and 21% when Hy/Hy6 and HcJ/HjS are 
used. 

6.2. Comparison between Balmer Decrements 

For a sample of 35 objects we have measurements of both 
the Hy/H/3 and the H6/H/3 ratios and can compare the at- 
tenuations derived using the two Balmer decrements (Figure 
|9|. The scatter between the two measurements for the same 
objects is significant. However, the mean values of the at- 
tenuation inferred separately from the Hy/H/? and Hd/H/S are 
consistent with the one-to-one relation, within the Icr mea- 
surement uncertainties (Figure |9]). The Pearson correlation 
coefficient for the sample is 0.33 and the Spearman correla- 
tion coefficient is 0.38. These correspond to 95% significance 
of the correlation. The mean (median) difference between the 
two measurements is 0.19 (0.33) +1.3 mag (1 cr), generally 
consistent with no systematic. 

6.3. Balmer Decrements in z — 0.8 Ha-selected Galaxies 

The depth of our IMACS spectroscopy enables the detec- 
tion of Hy6 and Hy lines with S/N sufficient for the determina- 
tion of the nebular attenuations for 16% of our spectroscopi- 
cally confirmed Ha emitters at z = 0.8 (N=55). The H/3 and 
H6 lines are detected with sufficient S/N for the analysis for 
15% of the sample (N=50). In this section, we first exam- 
ine the distribution of the observed Balmer line ratios of the 
z = 0.8 Ha emitters and the SDSS sample, and then com- 
pare the distributions of A(Ha) for the z - 0.8 and the SDSS 
samples. 

The observed Balmer line ratios for individual z - 0.8 
galaxies as a function of J-band magnitude, mass and 
Ha+[N11] flux are shown in Figure [TOl The Hy/HyS line ra- 
tios are plotted in the top row, and the H6fH/3 line ratios in 
the bottom row. The difference in the number of objects plot- 
ted in each panel is due to variations in the available data. The 
H6fH/3 sample is not entirely a subset of the Hy/H/3 sample. In 
some cases the Hy line is contaminated with sky emission but 
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Fig. 7. — Left: The observed J-band magnitude correlates tightly with the stellar mass of the galaxy for the Hff emitters in SXDS sample. This relation allows 
us to determine the masses of the SA-22 Ha emitters, for which PSF-matched broad band photometry has not yet been produced. Right: The MEx diagnostic 
diagram for the SA-22 sample, using the J magnitude - stellar mass relation allows us to identify three AGN (blue triangles). The local SDSS sample is also 
presented for comparison with symbols identical to those in Figure|6] 




Fig. 8. — Comparison between the attenuation A(H(!') derived from different Balmer line ratios in the SDSS sample as a function of J magnitude, stellar mass 
and SFR. Shown are the binned weighted means (open symbols, lines) and the weighted standard deviation in each bin. The bins are slightly olfset for clarity. 
The dashed grey horizontal line indicates A(Hq')=1 for reference. In the middle panel the range of masses of the NewHa sample is also indicated (vertical solid 
line). The Hafiifi and Hy/HyS ratios yi eld consist ent attenuation as a function of all parameters. The H5/H/3 ratio deviates at bright magnitudes, high masses and 
low star-formation rates (see Groves et al.l201 iT) . 



the W6 line is clear. Such objects are included in the H<5/Hj6 
sample but excluded from the Hy/H/? sample. 

Figure 10 shows a fair amount of scatter in the observed 
line ratios. First order linear fits to the data indicate relations 
in the expected directions: that the ratios decrease (i.e. atten- 
uation increases) towards brighter magnitudes, higher masses 
and higher SFR. However the data are equally consistent with 
a constant ratio as they are with a linear fit (based on the re- 
duced of the fits). The mean (median) line ratios in the total 
sample as a function of J-band magnitude (the largest sample) 
are Hy/H^^ 0.42(0.42) + 0.06 and mjVip^ 0.23(0.22) ± 0.06. 
When the AGN and composite galaxies are excluded, the 
mean line ratios and standard deviations are: 0.42 + 0.07 
and 0.22 + 0.04, respectively. Figure [10] also compares the 
z - 0.8 sample with the conditional distribution of the z = 
SDSS sample (i.e., each bin in the distribution is normalized 
separately). For the Hy/Hj6 sample, the distribution of the 
z - 0.8 line ratios is consistent with that of the z = sam- 
ple and the least squares fit to the z - 0.8 data agrees well 



with the 0.5 percentile local contour. The systematic shift in 
the SDSS H(5/Hj6 values towards higher Inie ratio (i.e., lower 
attenuation) as discussed in Section 6.1 creates an apparent 
difference between the SDSS weighted mean and the linear 
fit to the z - 0.8 data. The mean (median) line ratios (AGN 
included) correspond to an attenuation for the NewHa sam- 
ple of 0.8(0.9) + 1.0 (Hy/H^), 0.7(0.9) + 1.0 (H5/HyS) and 
0.8(0.9) +1.0 (joint) for the J-magnitude sample, which are 
consistent with the average attenuation from the Hy/HyS line 
ratio in the SDSS sample (1.0 ± 0.8). 

Some of the NewHo- objects exhibit line ratios above the 
Case B line which would indicate negative attenuation. Rather 
than having unphysical line ratios, these objects are likely 
scattered up due to errors in their line measurements. The off- 
sets above Case B are generally consistent with the measure- 
ment errors. The local SDSS sample also exhibits a spread of 
points which populate the area of the plot above Case B. 

The overall trends of A(Hq') based on the individual line 
measurements for the z - 0.8 and the z-0 samples are shown 
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Fig. 9. — Comparison between the values of the broadband attenuation A,, as 
derived from Hy/H/3 and from Hd/H/3 for a subsample of 34 objects where 
all three lines are detected at S/N> 5. Contours mark the locus of SDSS 
galaxies. The one-to-one line is shown (dotted line). The correlation between 
the two values is significant at the 95% level. 

in Figure [TTl again as a function of J-band magnitude, stellar 
mass and HQ'-f[Nll] flux. We show the binned weighted mean 
and weighted errors on the mean for the two samples where 
the measurement errors on the individual points are used as 
weightQ The NewHa points are binned such that each bin 
contains one third of the sample (blue points). The results 
from the Hy/H/S line ratio (top row) and the Hd/Hfi line ratio 
(bottom row) are both compared to the Hy/HyS-derived SDSS 
attenuation. Here we have excluded all NewHo- galaxies iden- 
tified as AGN and composite objects. In each panel we also 
show the total of the difference between the NewHo- and 
SDSS attenuation: 

2 _ y iAiHa,S PS S) - A(Ha,NewHa)f 
^ ~ ^ a-{A(Ha,SDSS)f + o-{A{Ha,NewHa)f 

In all panels the;^'^ is small indicating that the weighted means 
of the two samples are very similar. As in Figure [10] we also 
compare the reduced of a linear fit to the NewHtf binned 
points with the of a constant ratio. The linear fits in all 
panels have smaller than the constant ratio, indicating sig- 
nificant linear trends. These differences are smaller for the 
trend as a function of Ha-i-[NII] flux and larger for the trends 
as a function of J-magnitude and stellar mass. 

6.4. Balmer Decrements from Stacked Spectra 

In Section 6.2 Balmer decrement measurements for indi- 
vidual galaxies were examined as a function of J-band mag- 
nitude, stellar mass and Ha-H[NII] flux. However, the indi- 
vidual measurements exhibit a great deal of scatter, due to 
both intrinsic variations among the galaxies and to measure- 
ment errors. The measurement uncertainties stem from the 
limitations to proper continuum subtraction of the faint, high- 
order Balmer lines, and from measuring the line fluxes from 
low S/N spectra. The median error in the computed attenua- 
tion for individual galaxies in the Hy/HjS sample is 0.8 mag- 

We use the routine meanerr.pro which is adapted by Marc Buie from 
"Data Reduction and Error Analysis for the Physical Sciences", p. 76, Philip 
R. Bevington, McGraw Hill 



nitudes, and 0.7 magnitudes in the WSIWjS sample. The cor- 
responding errors in the SDSS sample are 0.6 and 0.3 magni- 
tudes for Hy/HjS and H(5/Hj6, respectively. However, whereas 
the SDSS shows statistically significant correlations with the 
three properties plotted (e.g., as discussed by IGarn & Besi 
|2OT0l) , the errors, combined with the relatively small size of 
the sample of z = 0.8 Ha emitters with sufficient S/N mea- 
surements of the Hj8, Hy, and H(J lines, obscures any potential 
trends. Furthermore, the objects with individual Balmer line 
measurements are likely a biased sample and will tend to have 
higher luminosities, higher SFRs, higher masses and/or lower 
attenuation. 

The last point is demonstrated in Figure [12] which com- 
pares the full sample of emitters with spectroscopic follow 
up and the sample of emitters with individual Balmer line 
detections. The objects with individual detections are more 
luminous, having brighter J-band magnitudes. The median 
J-band magnitude of the sample with individual line detec- 
tions i Jmed - 21.4) is 0.6 magnitudes brighter in J than the 
median for the full sample {J,„ed,fuii - 22.0). Faint galax- 
ies are underrepresented with only one object fainter than 
7 = 23 mag, with an individual line detection. The distribu- 
tion of Ha+[Nll]A65^4 fluxes is also skewed toward higher 
values: the median line flux for the full spectroscopic sam- 
ple is 13x10"'^ ergs s"' cm"^, while for the subsample with 
line detection the median is 27x10"'^ ergs s"' cm"^. As 
a result of contamination from sky lines at the positions of 
the Balmer lines, approximately half of the sample is ex- 
cluded from the stacks. However, the median J magnitude and 
Ha+[Nll]A65^4 flux of the remaining sample are very close 
to those of the full sample (J = 22.0 and //q'-h[NII]= 11.8) 
indicating that excluding the objects with sky line contami- 
nation does not bias the sample. Furthermore, a comparison 
between the distributions of the values of Ty (from the SED 
fits) for the included and excluded samples does not reveal a 
bias in the properties of the included sample. 

In order to address these issues, we create stacked spectra 
in bins of J-band magnitude, stellar mass and Ha-i-[NII]/}6584 
flux. The stacks contain both objects with and without indi- 
vidual detections of the Balmer lines. The bins are selected in 
a manner which achieves a minimum S/N of 10 in the li/3 line 
and a minimum of 10 objects per bin. This leads to uneven 
number of objects per bin (fewer objects in brighter bins) and 
uneven bins. The stacked spectra are created by co-adding the 
inverse-variance-weighted individual spectra. The error vec- 
tors are created by co-adding the individual inverse variance 
vectors. The final errors on the line measurements from the 
stacked spectra are the actual line measurement errors. The 
line measurements from the stacked spectra are compared to 
the error-weighted means of the local sample. The errors on 
the z = points represent the error-weighted error of the 
mean. 

The Balmer decrements measured from the spectra stacked 
as a function of magnitude are shown in Figure [T3] Out of 
the total sample of 341 objects at 0.78 < z < 0.83, the J- 
band stacks include 174 objects. AGN and AGN+SF objects 
are also exluded. The number of objects in each bin, the line 
fluxes, line S/N and attenuations in each bin are listed in Table 
[3] The points indicate a shallow decrease of the attenuation 
with increasing magnitude however the correlation does not 
have a high statistical significance - even without the brightest 
bin the Pearson coefficient is -0.31. The most luminous bin 
exhibits attenuation significantly below the local value. The 
large = 91 of the difference between the z=0.8 and z=0.0 
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Fig. 10. — Observed Balmer line ratios for individual z = 0.8 galaxies as a function of J-band magnitude (left column), stellar mass (middle column), and 
//q'+[NII]/16584 flux (right column). Objects with Hy/HyS line measurements are shown on the top row and those with Hi5/HyS line measurements on the bottom 
row. The numbers of star-forming galaxies (red symbols), galaxies with composite emission (cyan symbols) and galaxies hosting AGN (blue symbols) included 
in the analysis are indicated in the top left corners of the panels. The dotted lines indicate the ratios expected under the conditions of Case B recombination 
and the dashed lines indicate the least squares fit to the individual line ratios. In all cases, the slope of the least squares fit is in the expected direction - that the 
line ratios decrease with decreasing magnitude, increasing stellar mass and increasing //Q'+[NlI]/i6584 flux. Also shown is the conditional distribution for the 
local SDSS sample (0.05 to 0.95 percentile, shaded light gray). The local and z. = 0.8 samples occupy similar loci and the least squares fits (black, dashed line) 
agree well with the 0.5 percentile con tour (dark gray, solid line). The ofl'set of the local HS/Ufi sample towards higher line ratios is due to a systematic error and 
discussed in Section 6.1 and Groves et alJ i20Tl1) . 



samples is mostly driven by this bin. It is possible that the 
most luminous bin suffers from AGN contribution from ob- 
jects which cannot be individually identified as AGN because 
they have too low line S/N to be placed on the diagnostic dia- 
gram. The individual stacked spectra are shown in Figure [14] 
The Balmer lines of the stacked spectrum in the first bin are 
slightly broader than those of the fainter bins suggesting that 
AGN contamination may indeed be present. 

A number of studies in the local universe have shown 
that attenuation increases in star-forming galaxies with stel- 
lar mass, such that more massive galaxies are dustier. Figure 
[15] shows the stacked Balmer decrement measurements as a 
function of the stellar mass. A total of 168 individual spectra 
are combined in the stacks. The remaining objects are ex- 
cluded due to lack of mass measurements and prominent sky 
lines in the spectra. The attenuation shows a mild increase to- 
wards higher masses, more prominent in the Hy/H/S sample. 
Again, AGN contamination may be present in the most mas- 
sive bin where > 10.4, as expected from Figure [6] The 
= 7 is mainly dominated by this bin. The remaining mass 
bins closely follow the local relation. The stacked spectra are 
shown in Figure [16] and the line measurements for each bin 
are presented in Table [4] 

Finally, local studies indicate that dust con tent and atten- 
uation increase with the observed SFR (e.g., Hopkins et al.l 
l200lt[Sumvanetal.ll200T]:lGarn &Besll201 (J) so the same 
procedure is carried out as a function of Ha+[Nll]A65^4 flux. 
The results are presented in Figure [T7] 174 emitters are used 



in these stacks. The attenuation in the stacked spectra does 
slowly increase with //q'-h[NII]/16584 flux, in good agree- 
ment with the running weighted mean of the SDSS sample. 
The good agreement with the local sample is measured by the 
low = 4 of the difference between the two samples. The 
stacked spectra and line measurements are shown in Figure 
[18] and Table [5] respectively. 

Figures[T3][l5]and[T2]also show a comparison between the 
stacked spectra which include only star-forming galaxies and 
stacked spectra which also include the AGN and AGN-t-SF 
objects identified in Section 5. The line ratios of the stacked 
spectra which include AGN and composite galaxies are over- 
all consistent in most bins with the stacks which do not. How- 
ever, greater than Icr differences appear in the most massive 
and most luminous bins, such that the stacks which include 
AGN have lower attenuation. In the remaining bins AGN 
which were not individually identified should have a minor 
effect. 

7. DISCUSSION 

In the previous sections we presented Balmer decrement 
measurements for a sample of Ha-selected star-forming 
galaxies at z = 0.8. Both the individual and stacked spectra 
show that the average attenuation of the sample is A{Ha) ~ 1 
magnitude. The attenuation as a function of both stellar mass 
and SFR is similar to those at z = in comparison to the SDSS 
sample. In this section we compare these results to those from 
other local and intermediate-redshift samples and discuss the 
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Fig. 11. — Comparison between the weighted mean Balmer decrements of 7. = 0.8 star-forming galaxies and of the local SDSS sample contours as a function 
of J-band magnitude (left column), stellar mass (middle column), and //a+[NII]/l6584 flux (right column). Blue squares indicate the binned weighted mean and 
the weighted error of the mean for the 7, = 0.8 sample with the number of galaxies in each bin indicated. The individual measurements are also shown (small 
light blue symbols). For the local sample we also indicate the weighted mean and the standard error of the mean (dark gray line, light gray shading). The fiducial 
value of A(Hq') = 1 magnitude is also shown (dotted line). We display separately the attenuation derived from the Hy/H/? (top row) and from the H5/H/? (bottom 
row) decrements. The 7 = 0.8 sample exhibits trends consistent with those of the z = sample as indicated from the;('^ estimates in each panel. 

tency, despite the different methods for determining the atten- 
uation. 

The relation between SFR (or Ho- luminosity) and at- 
tenuation is frequently studied in light of efforts to derive 
attenuation-corrected SFR with minimu m auxiliary informa- 
tion. A num ber of previous studies (e.g., Hopkins et al. 200^ 
[SuUiyan eTal. 2001; Dopita et al. 2002; Perez-Gonzalez et al 
2003; Schmitt et al. 2006; Caputi et al. 2008) have found such 
a correlation locally, with the implication that galaxies with 
higher SFRs are also dustier Most such works have found a 
generally consistent slope and normalization. Here, we com- 
pare our results with the following two local relations: 



possible limitations of our analysis. 

7.1. Comparison with Other Extinction Relations 

Previous studies have found a correlat ion between mass 
and at tenuation in star-forming galaxies. iBrinchmann et al.l 
(120041) notice a clear increase in dust attenuation as a func- 
tion of mass as well as a wider range of attenuations for 
massive galaxies. For the SDSS sample they find that the 
mass-weighted attenuation A{Ha) is 1.3 with a spread of 
~ 0.7 magnitudes. Based also on SDSS and using a princi- 
pal component analysis, Garn & Best (2010) claim that mass 
is the most fundamental property which correlates with the 
attenuation, while SFR and metallicity are secondary, mainly 
brought about by the dependence of these parameters on mass. 
|Garn et al. (2010) compare SFR indicators for a sample of z = 
0.84 Ha-selected emitters similar to ours, estimating the at- 
tenuation from the ratio of the Ha and lAjim fluxes. They find 
that the mass-attenuation relation at that redshift is similar to 
the local one, with a possible plateau at M, > 1.6 x IO'^Mq. 
With the current sample we also find that z - 0.8 galaxies 
show a mass-attenuation relation similar to the local one. In 
Figure [Is] we comp are our results with the derived relation of 
iGarn & Best! (120101) . based on the SDSS sample: 



a- 
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logSFR = (0.597 + 0.021)^ + (-2.191 + 0.086), 
(Perez - Gonzalez et al. 2003, Eq.4) 



(4) 



SFRo^SFRi- 2.614 xlogn^ 



0.191 xSFRi + 3.^34] 



(Hopkins etal. 2001, Eq.5) 



(5) 



A(Ha) = 0.91 + O.llX + 0.1 IX^ - 0.09X' 



(3) 



where X = logj,)(M»/10'"Mo). We find agreement within the 
uncertainties, indicating little or no evolution in the amount 
of attenuation at a gi ven stellar mass ov er the last 6.8 Gyrs. A 
comparison with the Garn et al] (12010*) relation for their sam- 
ple of z = 0.84 emitters (their Figure 6b) also shows consis- 



The iPerez-Gonzalez et alj (12003b study is based on a sample 
of nearby objective-prism-selected star-forming objects. The 
SFR is derived from the observed Ha luminosity, using the 
iSullivan et al] (1200 Ih c onversion. The [H opkins et al. (200l]) 
relation is based on the IWang & Heckman (.1996) sample of 
z ~ FlR-detected galaxies and shown to be a good fit to a 
large (~ 500 objects) compilation of local star-forming galax- 
ies from iCram et ah (,1998.) . The relations between attenua- 
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Fig. 12. — Top row: Comparison between the full spectroscopic sample 
(dotted lines) and the Hy/HfS sub-sample (dashed lines) in terms of 7-band 
magnitude (first panel) and Ha+[T<lll]A65&4 flux (second panel). Also shown 
is the sample which excludes the spectra contaminated by sky lines (solid 
line). The objects in the HSjHfi sub-sample have brighter / magnitudes and 
larger //Q'+[Nll]/i6584 fluxes than the underlying sample. This means that 
the attenuation measurements based on this sub-sample will not be valid for 
low mass, lower-SFR galaxies, but preferentially biased to high mass, high- 
SFR galaxies. Bottom row: Same as the top row but for the HSjHfS sub- 
sample. The shift towards brighter / magnitudes and high //q'+[N11],16584 
fluxes is even more evident. 

tion and star-from ation have been extended up to z = 0.84 by 
iGarn et all (12010 ) (albeit using a different attenuation indica- 
tor), in agreement with local ones: 



Ana = -21.963 -H 0.562 logio 



-'Ha,ohs 



erg s 



-1 



(6) 



Ou r z = 0.8 sa mple is in agreement with both the local and 
the lGarn et al.l (12010) relation as shown in Figure [TT] 

Both comparisons show that the dust content of star- 
forming galaxies as a function of stellar mass and SFR have 
not significantly changed during the second half of the life- 
time of the Universe. This finding appears to be at odds with 
studies which have measured increased attenuation in galax- 
iesatO.5 <z< 1.0 (e .g, Flores et al. 2004; Liang et al. 2004; 
iRodrigues et al.l2008t:.Villar et al..2008;,Savaglio et al.,2005.) . 
However, the majority of such studies have targeted LIRGs 
and ULIRGs. While a detailed comparison between the at- 
tenuation of LIRG and non-LlRG galaxies within our sam- 
ple is beyond the scope of this paper, we can briefly exam- 
ine the issue here. Deep 24 fiva MlPS/Spitzer observations 
of a portion of our SXDS field have been acq uired by the 
Spitz er UKIDSS Ultra Deep Survey (SpUDS, iDunlop et alj 
120071) . reaching 3cr flux level of 25 mJy. This limit corre- 
sponds to a total IR luminosity (TIRL ) L]r = 2 x lO'** Lq 
for an object at z = 0.8 using the Riek e et al.l (l2009h relation 
between observed 24 //m flux and TIRL. A total of 50 Ha 
emitters with spectroscopic redshifts fall within the SpUDS 
footprint. We match 24 yt/m counterparts to 30 of these emit- 
ters (using a 3" matching radius) and find that 40% of the 
sample in the field (20 objects) have LlRG-like TIR luminosi- 
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Fig. 13. — Mean attenuation based on the Balmer decrements measured 
from stacked spectra. Stacked spectra are created in bins of apparent NEW- 
FIRM J-band magnitude. Each spectrum in the stack is weighted by the 
inverse variance. The top axis indicates the corresponding absolute J-band 
magnitude for the z = 0.8 sample. The stacking is performed for all ob- 
jects, including AGN and composite objects (open squares) and excluding 
both composite and AGN objects (filled red squares). All errors are based on 
the line measurement errors. The number of z = 0.8 spectra in each bin is 
indicated in red above the symbols. Also shown is the error-weighted mean 
attenuation for the local SDSS sample as a function of magnitude and the 
Icr weighted error on the mean (dark gray line, light gray shading). The 
horizontal dotted line indicates the fiducial value of A(Ha) = 1 magnitude. 
X~ = 91.42 is a measure of the differences between the two samples. The 
bins and the con'esponding hne fluxes are listed in Table|3] Figure n4l shows 
the binned spectra. 

ties (L/R > lO" Lq). Only five of the LIRGs have measured 
Balmer decrements so, instead of using the decrements, we 
derive A(Hq;) by comparing the SFR derived from the Ha flux 
and fr om the 24 yt/m flux in a manner similar to iGarn et all 
( I2OIOI) . We find that the median attenuation of the non-LIRG 
and LIRG galaxies is A(Ha)= 1.0 + 0.8 and 1.8 + 0.7 mag- 
nitudes, respectively, indicating a much higher attenuation for 
the LIRGs, consistent with prior studies. Within our sample, 
the LIRGs have brighter J-band magnitudes and dominate the 
sample for J < 21.5. The LIRG galaxies also have higher 
Hq'-h[N1I]/16584 fluxes: the median Ha-H[N11]/16584 flux for 
the LIRG sample is 2.7 times higher than the median of the 30 
24fj.m non-detections and non-LIRG sources: 20.5 x 10 vs. 
7.6 X lO^^^ergs s^^ cnr^. The bright J-band magnitudes and 
elevated median Ha flux indicate that the LIRGs in our sam- 
ple occupy the parameter space where we have already found 
elevated attenuation (Figures [T3] and [TsTl. and have a smaller 
contribution to our sample at low-luminosity and low-SFR. 

However, elevated attenuations are also found in studies of 
galaxies with lower SFRs. For example, in a study of 69 
0.4 < z < 1. galaxies in the spec tr oscopic Gemin i Deep 
Deep Survey dAbraham et al.1 120041) . ISavagho et al] (120051) 
findHy/Hy3= 0.29 + 0.03 andH(5/Hj8= 0.18 + 0.02 in a stacked 
spectrum of all objects, lower than most individual detec- 
tions in our sample (corresponding to A(Hq') of 3.34 ± 0.2 
and 1.68 ± 0.15, respectively). Villaretal. (2008) also find 
higher attenuation in their narrow-band Ha-selected z - 0.84 
sample - A(Ha) - 1 .48. While our results are consistent with 
this value within the errors, we do not find a shift in the aver- 
age attenuation of galaxies with physical properties similar to 
the local sample. The difference may be due to the different 
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Fig. 14. — Weighted stacked spectra, divided in six bins by J-band apparent magnitude. Shown are the stacked spectra (gray line), the continuum fits (red 
dotted line) and the continuum subtracted spectrum (offset by 0.5 vertically, black line). The magnitude increases top to bottom. Vertical dotted lines indicate the 
positions of the Balmer fines US (left), Hy (center) and Hfl (right). 
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Fig. 15. — Mean attenuation based on the Balmer decrements measured from 
stacked spectra, binned by the logarithm of the stellar mass. The symbols are 
the same as in Figure [731 The local relation between attenuation and stellar 
mass I Ga rn & B est 2010) is also plotted (dot-dash line). The weighted mean 
attenuation for the local Hy/HyS SDSS sample as a function of stellar mass is 
shown (solid dark gray line) with errors indicating the weighted error on the 
mean (Ught gray shading), = 6.79 is a measure of the differences between 
the two samples. Figure ]l6l shows the binned spectra and Table|5]lists the Une 
fluxes and attenuation values for each bin. 

attenuation measurement method - instead of Balmer decre- 
ments Villar et al. (2008) use determine t he mean at tenuation 
for their sample via the ¥dustl^Fuv ratio (Buat et al. 2005). If 
the FUV continuum, the Balmer lines and the dust emission 
sample different physical scales, the attenuations may differ 
Nevertheless, a comparison is valuable, because their sample 
selection is very similar to ours. Similar to their work, we 
derive a single attenuation for our full sample, to determine if 
their average is weighted towards more massive and luminous 
(and more attenuated) sources. However, even if, instead of 
splitting our sample in bins by SFR, we combine them in a 
single stacked spectrum, we find lower, not higher values of 
the attenuation: A(Ha)=0.25 + 0.75 from the ratio of Hy to 
H/3 and A(Ha')= 0.55+0.65 from the ratio of to Hy3. These 
values are ~ 0.5 magnitude lower than the median attenuation 
from the individual detections. They are also > 1 magnitude 
lower than the values of Villar et al. (2008) and inconsisten t 
within 1 cr with the measurements by Savag ho et all (l2()05h . 
The stacked spectrum is dominated by low-luminosity galax- 
ies which in both SDSS and in NewHa have attenuation of 
< 1 magnitu de. Even if we wei gh equally all galaxies in 
the manner of ISavaglio et al.l (12005 ) who normalize all spec- 
tra to one before stacking them, we still measure attenuation 
of A(Ha)~ 1.0 magnitude, inconsistent with their results. In 
conclusion, we do not find evidence for systematically ele- 
vated attenuation relative to the local sample. 

7.2. Limitations of the Balmer Decrement Analysis 

Problems may exist with our approach to determining the 
attenuation using the high order Balmer decrements. In this 
section we discuss the limitations of using Hy6, Hy and W6 to 
constrain the attenuation. 

The high-order Balmer lines are faint. Hy is 5.85 times 
fainter than Ha (assuming Case B), and WS another factor of 
two fainter We have shown that this can be overcome by deep 
spectroscopic observations, relatively high spectral resolution 



and stacking of sources. Assuming Case B for the Hy/HyS ra- 
tio, we can calculate that in order to achieve an error in A(HQf) 
of 0.5 magnitudes both lines must be detected at a S/N> 20. 
To achieve the same precision with H(5/Hy6 we need S /N > 14. 
In comparison, the median error of the individual detections 
of our sample is 0.7 magnitudes. There are only seven ob- 
jects with Hy detected with S/N> 20 and six objects with ¥16 
detected with S/N> 14, all of them identified as AGN or com- 
posite objects in Section 5. The mean Balmer ratios for these 
two samples ai-e 0.42 + 0.04 and 0.22 + 0.03 for Hy/HyS and 
H5/Hy6, respectively, consistent with our findings for the full 
samples. Furthermore, by stacking individual spectra we have 
been able to measure the lines with S/N> 15, again achieving 
consistent results. 

A complementary line of work would be to construct the 
Ha/Hy6 line ratio using the Ha fluxes from our NEWFIRM 
narrow -band NB118 photometry. Assuming that the contin- 
uum noise is the same at all wavelengths, the S/N in the Ha 
line will be ~ 14 times higher than that of H/5 and ~ 30 
times higher than that of Hy (Case B). Such an analysis re- 
quires a careful examination of the possible sources of sys- 
tematic error when combining spectroscopic and photometric 
data such as differences between the apertures within which 
the H/3 and Ha line fluxes are measured. Furthermore, the 
NBl 18 bandpass is wide enough to include flux from the [N 
II] /i/i6548, 6583 emission lines for the Ha emitters and we 
need to correct the Ha fluxes for the contamination in or- 
der to construct the Ha/H/3 decrement. The Ha/[NII]/16584 
ratio is a metallicity indicator and thus varies as a function 
of equivalent width (i Villar et al. 2 008). B-band luminosity 
dKennicutt et al.l2008l:lLee et al.l200 9) and mass, but the scat- 
ter in these relations is significant. Near-IR spectroscopy is 
required to directly measure the Ha/[NII]/16584 ratio. 

In addition to the measurement sources of error, there are 
physical weaknesses in using Balmer decrements to deter- 
mine attenuation. One regime where the optical attenuation 
measurements break down is in extremely dusty galaxies. 
Even though many such objects may exhibit an emission line 
spectrum indicating moderate attenuation, the nebular redden- 
ing may bear little resemblance to the actual dust content. In 
such cases the emission comes from the least-obscured re- 
gions of the galaxy and is not representative of the galaxy as 
a whole. This issue is manifested in observed discrepancy 
between the Balmer decrement attenuation and IR-to-Ha at- 
tenuation. In such galaxies, the Balmer decrements underesti- 
mate the actual attenuation. The mid-IR and far-IR emission 
is much more sensitive in unmasking obscured star formation 
and can be used to address this issue in future work. 

8. SUMMARY AND FUTURE WORK 

We have examined the dust attenuation in normal star form- 
ing galaxies at z = 0.8 using the Hj/H/S and H6/H/3 line ratios. 
We have carried out deep spectroscopic observations for a 
sample of narrow -band Ha selected galaxies, which allows us 
to detect H/3, Hy and H6 with S/N> 5. We apply several dif- 
ferent techniques to isolate AGN and composite objects from 
the sample: the MEx diagnostic, the UV variability and the 
spectral type. In total, we identify 17 AGN and seven com- 
posite objects which indicates a strict lower limit on the AGN 
fraction of 5%, in line with other studies. 

We find good agreement between the distribution of atten- 
uations in local galaxies and those from the NewHa sample. 
The mean attenuation of galaxies with individual line detec- 
tions is A(Ha) = 0.8 + 1.0 and 0.7 + 1.0 for the Hy/HyS and 
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Fig. 16. — Same as Figure [T4l but for the logarithmic mass bins. 
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in the spectra and probe a more representative sample of the 
overall star-forming population than possible with individual 
detections. The stacked spectra also indicate consistency be- 
tween the attenuation at z = 0.8 and z = for galaxies with 
the same mass and SFR. 

The analysis presented here is solely focused on the Balmer 
decrement as an attenuation indicator. As we have pointed 
out, this method is not applicable for very dusty objects where 
the nebular emission originates only in the least dusty, op- 
tically thin regions of the galaxy. In order to explore this 
regime, we plan to use the MIPS 24/im observations from the 
SpUDS survey which overlaps with part of our SXDS field. 
Further work will also include comparison of the Balmer 
decrement results with measurements of attenuation from the 
UV slope and SED models. 
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APPENDIX 

APPENDIX: BALMER DECREMENTS FROM OPTICAL SPECTROSCOPY 

The Balmer emission lines arise from the recombination and subsequent cascade of electrons to the n = 2 level of the hydrogen 
atoms in the interstellar medium. The preceding photoionization is caused by the radiation field of hot, young, massive stars. 
Collisional excitation does occur as well (e.g., Ferland et al. 2009) but photoionization and recombination are the dominant 
processes. 

The simple atomic structure of hydrogen has made it possible to determine the expected rates of transition and the emission 
line ratios as a function of the temperature and electron density of the interstellar medium. Two general cases are considered 
when discussing the Balmer line emissons. In the Case A recombination, the nebula is optically thin to the photons of the Lyman 
series (i.e., photons emitted from transitions to the « = 1 level of hydrogen). Alternatively, in the Case B recombination model, 
the nebula is optically thick to photons of the Lyman series but optically thin to Lya and photons of the higher hydrogen series. 
Case B is typically assumed for determining the intrinsic line ratio. Here we assume Case B recombination and electron density 
and temperature typic al for HII regions: 10"* cm"^ and 10,000 K. The expected intrinsic Balmer line ratios for these condition 
are (lOsterbrock 111981 : 

(Ha/Hj3)i„,^2M (Al) 
(Hy/H/3)i,„ ^0.469 (A2) 
(H6/H/3)i,„ ^0.260 (A3) 

The ratios are only weakly sensitive to electron density and more sensitive to temperature. For example the Ha/H/3 ratio varies 
by only 0.05 (2.86 to 2.81) f or electron densities «c = 10^ to lO*" cm'^ but by 0.29 for temperatures between 5000K and 20,000K 
dPopita & Sutherlandl2003h . However these variations are still small relative to the effects of dust. 

Our optical spectra only span the 3500 to 5300A region and do not contain the Ha emission line at the redshift of the sample. 
We do have Ha observations from our NBl 190 photometry but the direct comparison between the spectroscopic and photometric 
line fluxes is impeded by the need of aperture corrections which we defer to future work. In this paper we focus on the Hy/H/3 
and H6/H/3 decrements. The effect of dust is to diminish the observed emission of the short-wavelengths more than the long 
wavelengths. The attenuated observed line ratios will then be smaller than the intrinsic ones due to the effect of dust. The color 
excess due to this dust attenuation can be expressed as: 
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Fig. 17. — Mean attenuation based on the Balmer decrements measured from 
stacked spectra, bi nne d by their //q'+[N1I]/16584 flux. The symbols are the 
same as in Figure \\3\ The weighted mean attenuation for the local Hy/HyS 
SDSS sample as a function of stellar mass is shown (solid dark gray line) 
with errors indicating the weighted error on the mean (light gray shading). 
X'^ = 4.31 is a measure of the dift'eren c es bet ween the two samples. The 
local relations of [Perez- Gonzal ez et all {2003h . [Hopkins et al. ( 2001), and 
IGam et alj tICTTOl) as well as the z = relation of Gam & Best (2010) are 
also included as indicated by the legend. Figure [TSl shows the binned spectra 
and Table[5]lists the line fluxes and attenuation values for each bin. 

H(5/HjS samples respectively. The mean (median) of the joint 
sample is A(Hq')= 0.8(0.9) + 1.0. 

We stack the spectra in bins of J-band magnitude, stellar 
mass and Hq;-i-[NII]/16584 flux in order to increase the S/N 
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Fig. 18. — Same as Figure [T4l but for the binning in NBl 190 flux in units of 10 " erg s ' cm ^. 
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A{Hli)-A{Hy)^ 
= -//r) = -2.5 xlogio 



{HylHpUs 



(A4) 
(A5) 



where {Hy/Hf})ohs denotes the observed Hne ratio. This color excess can be related to the broadband color excess E{B — V) via 
an attenuation curve: 



A(A) = k(A)E(B - V), 



(A6) 



where k( A) is the value of the attenuation curve at the wavelength A. Using the attenuation curve we can express the attenuation 
at Hy and H/3 in the equation above: 



E(H/3 - Hy) = E(B - V)[k(H/3) - K(Hy)] 



(A7) 



The functional form of the attenuation curve is a matter of choice. We follow the common practice and adopt the lCalze"tt]| 
( 119971) attenuation curve in the form: 



/(i) = 2.659(-1.857+ 1.040jt:) + «;„ 
for 0.63//m < A < 2.20yum; 
= 2.659(-2.156 + 1.509x - 0.198x2 + 0.01 Ix^) + R'„ 
for 0.12jUOT < A < 0.63yum, 



(A8) 
(A9) 
(AlO) 
(All) 
(A12) 



where x - l/Ais the wave number and R[, - A(V)/E(B - V) - 4.05 is the ratio of the total to selective attenuation. Using this 
prescriptions, we find K{Ha) - 3.33, k(H/3) - 4.60 and K(Hy) - 5.12 and k(H6) - 5.39. The expression for E{H/3 - Hy) in 
Equation 7 can be then substituted in Equation 5, reaching a solution for the broadband color excess E(B - V): 



EiB - y) = 



E{H/3 - Hy) 
k(H/3) - K(Hy) 

-2.5 
KiH/5) - K{Hy) 



X log 10 



0.469 



{Hy/H/3) 



obs 



Similarly, the broadband color excess can be expressed as a function of the H6/H/3 Balmer decrement is: 



(A13) 
(A14) 



E(B - V) ■■ 



E(H/3 - H6) 
KiHp) - k(H5) 

-2.5 
k{HI3) - k(H6) 



X log 10 



0.260 



(H6/Hf5)ob. 



(A15) 



(A16) 



From here we can express the attenuation at an arbitrary wavelength via Eq. IA6I Most commonly the attenuation is determined at 
Ay = 5500A, where k(V) - 4.05. For the purposes of our work, we are also determine A(Ha), the attenuation at Ha {A = 6563A), 
where /((//a) = 3.33. 
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Nebular Attenuation at z = 0.8 



TABLE 2 

AGN AND Composite Candidates 



NewHfj- 
Field 




RA 


Dec 


FUV 
Variability 


Spectral 
Type 


MEx" 


Comment 


SXDSN 


02 


18:22.5 


-04:30:36.1 


+ 


+ 


- 


- 


SXDSS 


02 


17:42.6 


-05:03:51.4 


- 


+ 


- 


BL* 


SXDSW 


02 


17:01.3 


-04:54:52.6 


- 


+ 


- 


- 


SXDSW 


02 


16:14.7 


-04:46:05.3 


- 


+ 


- 


BL AGN 


SXDSN 


02 


17:43.0 


-04:36:25.0 


- 


+ 


- 


BL AGN 


SXDSS 


02 


17:21.9 


-05:03:49.2 


- 


- 


Sy2 (0.412)/L(0.324) 


- 


SXDSS 


02 


16:52.5 


-05:22:30.5 


- 


- 


L (0.487) 


- 


SXDSS 


02 


18:20.1 


-05:11:53.3 


- 


- 


L(0.5)/Sy2 (0.5) 


- 


SXDSS 


02 


16:50.9 


-05:21:15.1 


- 


- 


L (0.561) 


- 


SXDSW 


02 


17:09.9 


-04:55:56.7 


- 


- 


Sy2 (0.667) 


- 


SXDSW 


02 


16:31.7 


-04:57:25.5 


- 


- 


Sy2 (0.857) 


- 


SXDSW 


02 


15:51.8 


-05:09:47.7 


- 


- 


L(0.341)/Sy2 (0.317) 


- 


SXDSW 


02 


17:16.3 


-04:59:57.4 






L (0.519) 




SXDSN 


02 


18:16.6 


-04:26:34.4 






Comp (0.5 )/Sy2 (0.5) 




SXDSS 


02 


17:37.1 


-05:06:22.2 




+ 


Comp (0.610) 




SXDSS 


02 


17:21.8 


-05:17:24.4 






Comp (0.504) 




SXDSS 


02 


16:45.2 


-05:18:03.5 




+ 


Comp (0.570) 




SXDSS 


-05:16:48.4 


-05:18:42.5 






Comp (0.560) 




SXDSW 


02 


16:36.8 


-05:05:54.0 






Comp (0.802) 




SXDSW 


02 


15:48.4 


-04:53:55.9 






Comp (0.723) 




SXDSW 


02 


16:57.9 


-05:05:39.5 






Comp (0.467) 




SA22 


22 


17:30.8 


+00:12:18.9 






+ 




SA22 


22 


17:32.5 


+00:21:33.4 






+ 




SA22 


22 


17:12.0 


+00:12:44.1 






+ 





^ Classification according to 'Juneau et ali 1201 11, see Figure |6). The designations indicate LINER (L), 
Seyfert Type 2 (Sy2) and Composite SF+AGN (Comp). The fractional probability that the object belongs 
to a given class is given in brackets. Objects are assigned the class which has the largest fractional 
probability. The probabilities of the four classes (L, Sy2. Comp and SF) add up to 1. 
^ Broad line. 
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TABLE 3 

Stacked Spectra as a Function of J-band Magnitude ^ 



Bin J[mag] N'' Mean S/N'^ U/S^ AH/3 Hy AHy H<5 AUS A(Hy/Efi) A(Hy/H8) 



1 19.45 -20.75 16 3.7 3.12 0.15 1.65 0.12 0.95 0.09 -0.82 -1.1 

2 20.75-21.25 24 2.8 3.31 0.15 1.29 0.12 0.71 0.09 1.25 1.29 

3 21.25-21.75 28 2.1 3.45 0.16 1.57 0.10 0.78 0.08 0.20 0.93 

4 21.75-22.25 33 1.6 2.33 0.14 1.03 0.10 0.52 0.07 0.43 0.99 

5 22.25 -22.80 35 1.2 2.19 0.13 0.94 0.10 0.44 0.08 0.60 1.82 

6 22.80 - 25.00 38 0.9 2.00 0.12 0.84 0.10 0.48 0.07 0.75 0.53 



" The data presented here are for the stacked spectra which exclude individually-identified AGN and composite sources. 

Number of individual spectra used to create the stacked spectrum. 
" Mean continuum signal-to-noise in the stacked spectrum. 

^ Line fluxes are in units of 10"'^ ergs s"' cm"^. They are not corrected for attenuation or aperture losses. 



TABLE 4 

Stacked Spectra as a Function of Stellar Mass^ 



Bin M. N* MeanS/N'= Ufi^ AUfl Hy AHy H<5 AH<5 A(Hy/HS) A(Hy/HS) 



1 


8.00 - 9.25 


18 


0.8 


1.48 


0.14 


0.69 


0.10 


0.37 


0.08 


0.10 


0.22 


2 


9.25 - 9.75 


74 


1.2 


2.01 


0.13 


0.94 


0.10 


0.48 


0.07 


0.03 


0.62 


3 


9.75 - 10.10 


34 


1.9 


3.30 


0.16 


1.30 


0.10 


0.70 


0.08 


1.19 


1.41 


4 


10.10 - 10.40 


23 


2.4 


2.86 


0.16 


1.10 


0.11 


0.68 


0.09 


1.39 


0.63 


5 


10.40-11.15 


19 


3.3 


2.82 


0.16 


1.39 


0.13 


0.68 


0.10 


-0.35 


0.53 



^ The data presented here are for the stacked spectra which exclude individuaUy-identrfied AGN and composite sources. 

Number of individual spectra used to create the stacked spectrum. 

Mean continuum signal-to-noise in the stacked spectrum. 
^ Line fluxes are in units of 10"'^ ergs s"' cm"^. They are not corrected for attenuation or aperture losses. 



TABLES 

Stacked Spectra as a Function of Ha+[Nll]A65S4 Flux 



Bin Ha-H[Nn]'' Mean S/W* H/}"" AH/3 Hy AHy US AH<5 A(Hy/H^) A(Hy/H^) 



1 1-6 29 1.0 1.50 0.12 0.63 0.10 0.41 0.08 0.83 -0.29 

2 6 - 8 24 1.1 1.62 0.14 0.74 0.10 0.42 0.07 0.16 -0.01 

3 8- 12 31 1.34 2.66 0.13 1.07 0.10 0.49 0.08 1.08 2.31 

4 12-20 49 2.0 2.60 0.14 1.05 0.10 0.57 0.08 1.01 1.14 

5 20-40 25 2.5 3.68 0.15 1.61 0.11 0.86 0.08 0.47 0.78 

6 40- 120 16 3.7 5.04 0.15 2.22 0.10 1.00 0.08 0.41 1.88 



^ The data presented here are for the stacked spectra which exclude individually-identifled AGN and composite 
sources. 

* Line fluxes are in units of 10"'^ ergs s~' cm~^. They are not corrected for attenuation or aperture losses. 
Number of individual spectra used to create the stacked spectrum. 
Mean continuum signal-to-noise in the stacked spectrum. 



